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Abstract 

A  conceptual  transform  domain  communication  system  (TDCS)  is  shown  ca¬ 
pable  of  operating  successfully  using  M-Ary  phase  shift  keying  (MPSK)  data  mod¬ 
ulation  in  a  multiple  access  environment.  Using  spectral  encoding,  the  conceptual 
TDCS  provides  an  effective  means  for  mitigating  interference  affects  while  achiev¬ 
ing  multiple  access  communications.  The  use  of  transform  domain  processing  with 
MPSK  data  modulation  (TD-MPSK)  provides  higher  spectral  efficiency  relative  to 
other  modulation  techniques  (antipodal  signaling  and  cyclic  shift  keying)  considered 
previously  for  TDCS  applications.  The  proposed  TD-MPSK  technique  uses  spec¬ 
tral  encoding  for  both  data  and  multiple  access  phase  modulations.  Demodulation 
of  the  spectrally  encoded  TD-MPSK  communication  symbols  is  accomplished  using 
conventional,  multi-channel  time  domain  correlation  techniques. 

Analytic  expressions  for  TD-MPSK  probability  of  symbol  error  (Pe)  and  prob¬ 
ability  of  bit  error  ( Pb )  are  derived  and  validated  using  simulated  results  over  the 
range  of  signal-to-noise  ratios  typically  considered  for  communications.  This  valida¬ 
tion  includes  scenarios  with:  1)  multiple  access  interference,  2)  spectral  notching,  3) 
jamming  present  and  4)  combinations  of  all  three.  For  a  J/S  of  3.14  dB  and  a  E^/Nq 
of  6  dB,  Pb  dropped  by  up  to  a  factor  of  3  for  TD-QPSK  in  a  MA  environment  for 
the  case  when  spectral  notching  was  present  versus  the  case  when  spectral  notching 
wasn’t  present.  The  cross-correlation  between  communication  symbols  of  different 
synchronous  users  can  be  made  identically  zero  through  proper  selection  of  multiple 
access  phase  codes  (orthogonal  signaling).  For  a  synchronous  network  containing 
orthogonal  users,  Pe  and  Pb  are  unaffected  as  the  number  of  orthogonal  network 
users  increases.  For  a  J/S  of  3.14  dB  and  a  Eb/N0  of  6  dB,  Pb  dropped  by  a  factor 
of  12  for  TD-QPSK  in  a  MA  environment  for  the  case  when  spectral  notching  was 
present  versus  the  case  when  spectral  notching  wasn’t  present. 
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INTERFERENCE  SUPPRESSION  IN  MULTIPLE  ACCESS 
COMMUNICATIONS  USING  M-ARY  PHASE  SHIFT  KEYING 
GENERATED  VIA  SPECTRAL  ENCODING 

1.  Introduction 

1.1  Introduction 

Digital  communication  systems  have  become  an  integral  part  of  commercial 
day-to-day  operations  and  are  extensively  employed  across  today’s  modern  battle¬ 
field.  Military  applications  include  airborne,  seaborne,  and  ground-based  commu¬ 
nications,  navigation  such  as  the  Global  Positioning  System  (GPS),  and  unmanned 
aerial  vehicle  (UAV)  communications  and  control.  However,  because  of  the  nature  of 
war  itself,  military  applications  often  require  special  capabilities  not  commonly  asso¬ 
ciated  with  commercial  applications  communications,  e.g.,  secure  operation,  hostile 
(intentional)  radio  frequency  (RF)  interference,  and  transmission  of  large  amounts  of 
data  between  systems  (platforms)  which  often  times  were  designed  and  built  without 
any  need  or  forethought  of  integration. 

There  are  two  aspects  to  providing  secure  communications.  The  first  aspect 
involves  preventing  the  enemy  from  intercepting  the  transmitted  signal  (detecting 
that  communication  is  occurring).  The  second  aspect  involves  preventing  the  enemy 
from  deciphering  the  transmitted  signal  if  it  is  intercepted  (extracting  the  informa¬ 
tion  from  the  communications). 

There  are  several  types  of  RF  interference  that  occur  on  the  battlefield.  The 
first  type  of  interference  is  the  result  of  thermal  noise  and  is  inherent  in  all  com¬ 
munication  systems,  including  commercial  ones.  The  random  electron  motion  which 
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occurs  in  all  electronic  devices  induces  thermal  noise  [1].  The  second  type  of  bat¬ 
tlefield  interference  comes  from  additional  communication  transmitters  operating  in 
the  environment,  including  radio  stations,  television  stations,  and  cell  phone  sys¬ 
tems.  The  third  type  of  battlefield  interference  is  intentional,  i.e.,  communication 
and  radar  jammers  which  are  introduced  solely  for  the  purpose  of  disrupting  effective 
communications. 

The  transmission  of  large  amounts  of  data  allows  accurate  instructions  to  be 
transmitted  on  the  battlefield.  One  example  is  the  control  of  a  UAV  which  requires 
precise  instructions  from  its  home  base.  Another  example  is  a  smart  bomb  which 
requires  precise  targeting  instructions  from  an  aircraft.  A  reduction  in  transmission 
rate  from  the  communication  system  could  hinder  the  ability  of  these  weapon  systems 
to  effectively  accomplish  their  mission. 

In  communication  systems,  a  transmitter  sends  a  signal  over  a  channel.  By  the 
time  the  signal  reaches  the  receiver,  the  received  signal  consists  of  the  transmitted 
signal  plus  interference.  Traditional  communication  systems  deal  with  channel  in¬ 
terference  in  one  of  two  ways.  First,  the  receiver  may  notch  out  frequencies  where  it 
senses  interference  is  occurring.  Although  this  receiver  notching  effectively  removes 
the  source  of  interference,  it  also  removes  some  of  the  desired  signal  energy  which 
was  transmitted  in  the  signal.  Second,  the  receiver  can  estimate  the  interference 
and  subtract  the  estimate  from  the  received  signal.  The  main  problem  with  this 
interference  estimation  method  is  that  portions  of  the  transmitted  signal  could  be 
confused  with  interference  and  be  subtracted  from  the  received  signal  as  well  [2] .  In 
both  of  these  receiver  processing  techniques,  desired  signal  energy  is  removed  prior 
to  detection  and  estimation. 

Communication  systems  can  provide  secure  communications  by  encrypting  in¬ 
formation.  One  problem  with  encrypting  information  is  that  system  bandwidth  typi¬ 
cally  increases  (assuming  the  desired  data  rate  of  the  information  remains  the  same). 
Another  problem  is  that  the  signal  may  still  be  non-cooperatively  intercepted  even 
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though  it  may  be  undecipherable.  Although  the  information  content  of  messages  are 
not  discernable,  the  detection  of  a  signal  betrays  the  presence  of  the  communication 
system  and  could  still  lead  to  system  geolocation. 

Traditional  communication  systems  provide  high  data  transmission  rates  by 
efficiently  using  bandwidth  [3].  Different  data  modulation  techniques  permit  differ¬ 
ent  transmission  rates  and  provide  different  probability  of  bit  error.  Choosing  the 
appropriate  data  modulation  technique  generally  requires  a  trading-off  analysis  of 
data  transmission  rate  and  bit  error  probability. 

Spread  spectrum  systems  were  primarily  developed  by  the  military  as  a  method 
for  providing  secure  communications  in  environments  containing  significant  amounts 
of  RF  interference.  Two  widely  used  spread  spectrum  modulation  schemes  are 
direct  sequence  spread  spectrum  (DSSS)  and  frequency  hopping  spread  spectrum 
(FHSS)  [4],  FHSS  systems  mitigate  interference  effects  by  spectrally  hopping  from 
one  frequency  band  to  another.  If  interference  is  present  during  operations,  it  only 
affects  the  received  signal  for  the  portion  of  time  when  the  hop  frequency  is  in  the 
same  frequency  band  as  the  interference.  Frequency  hopping  systems  provide  secure 
communications  since  the  receiver  needs  to  know  the  frequency  hopping  sequence  to 
demodulate  the  received  signal  [4], 

In  DSSS  communications,  interference  is  suppressed  by  having  the  transmitter 
spread  the  desired  signal  over  a  bandwidth  which  is  much  wider  than  otherwise 
required.  The  receiver  then  despreads  the  received  signal.  Any  interference  present 
during  despreading  has  less  of  an  effect  on  received  signal-to-noise  ratio  since  it  spans 
a  smaller  fraction  of  the  received  signal  versus  a  traditional  communication  system. 
The  DSSS  technique  also  provides  security  since  the  intended  receiver  needs  to  know 
the  despreading  code  to  effectively  demodulate  the  received  signal  [4], 

The  transform  domain  communication  system  (TDCS)  uses  a  form  of  spec¬ 
tral  coding  to  take  interference  suppression  (avoidance)  one  step  beyond  traditional 
spread  spectrum  techniques.  The  TDCS  spreads  the  transmitted  signal  over  a  large 
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bandwidth  while  excluding  spectral  regions  containing  significant  amounts  of  inter¬ 
ference.  The  receiver  then  receives  the  transmitted  TDCS  signal  energy  without 
significant  amounts  of  interference  [5].  Because  of  the  spectral  coding  method  used 
by  the  TDCS  to  spread  the  signal,  the  time  domain  signal  appears  noise-like  and 
unintended  receivers  have  a  more  difficult  time  intercepting  it.  If  the  TDCS  signal  is 
non-cooperatively  intercepted,  the  signal  remains  nearly  impossible  to  demodulate 
without  prior  knowledge  of  the  spreading  code  used  to  modulate  the  signal. 

The  origin  of  the  TDCS  concept  can  be  traced  back  to  work  by  German  and  An- 
dren.  German  proposed  designing  a  direct  sequence  system  which  avoided  jammed 
frequencies.  This  was  done  by  having  the  transmitter  generate  a  spreading  code 
which  did  not  have  any  frequency  components  in  areas  containing  spectral  interfer¬ 
ence.  Andren  subsequently  designed  a  system  where  the  transmitter  generated  a  sig¬ 
nal  which  avoided  frequency  components  in  spectral  regions  containing  interference. 
The  receiver  only  extracts  energy  from  spectral  regions  containing  no  interference 
(i.e.,  using  only  spectral  components  where  the  transmitter  placed  energy)  and  then 
processes  the  signal  [5]. 

Radcliffe  combined  the  ideas  of  German  and  Andren  to  generate  the  conceptual 
TDCS  considered  here.  In  Radcliffc’s  TDCS  design,  time  domain  “basis  functions” 
(called  fundamental  signaling  waveforms  in  this  work)  are  generated  which  don’t 
contain  energy  in  spectral  regions  containing  interference.  The  TDCS  technique  then 
data  modulates  these  basis  functions  to  produce  the  digitally  encoded  waveform.  The 
types  of  data  modulation  considered  by  Radcliffe  were  cyclic  shift  keying  (CSK)  and 
binary  phase  shift  keying  modulation  (antipodal  signaling)  [5]. 

Swackhammer  expanded  Radcliffe’s  work  by  investigating  a  multiple  access 
version  of  the  TDCS.  Multiple  access  capability  allows  multiple  transmitter /receiver 
pairs  to  simultaneously  communicate  using  the  same  portion  of  the  RF  spectrum. 
Swackhammer  implemented  a  multiple  access  design  using  spreading  codes  (multiple 
access  codes)  generated  from  linear  feedback  shift  registers.  The  spreading  codes 
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were  used  in  a  phase-mapping  process  to  generate  a  set  of  unique  phase  values  to 
spectrally  encode  multiple  access  users.  The  types  of  data  modulation  considered 
for  his  design  were  CSK  and  BPSK  modulation  [6]. 

1.2  Problem  Statement 

Previous  research  has  shown  that  the  TDCS  technique  can  be  used  to  pro¬ 
vide  reliable  communications  in  hostile  RF  environments  while  generating  a  low 
probability  of  intercept  signal.  Previous  research  has  also  shown  that  multiple  ac¬ 
cess  capability  is  achievable  through  properly  phase  coding  (spectrally  encoding)  the 
spectrally  notched,  interference  avoiding  waveforms.  Previous  TDCS  research  pre¬ 
dominantly  focused  on  CSK  and  BPSK  data  modulation  for  communicating.  This 
work  introduces  a  form  of  spectral  encoding  which  provides  M-Ary  phase  shift  key¬ 
ing  data  modulation  and  orthogonal  signaling  for  multiple  access  capability.  As  with 
traditional  MPSK  signaling,  the  spectrally  encoded  waveforms  generated  here  can 
carry  more  information  than  BPSK  waveforms  using  the  same  bandwidth.  Perfor¬ 
mance  of  the  proposed  TD-MPSK  technique  will  be  first  characterized  for  a  single 
transmitter-receiver  pair  operating  in  an  environment  containing  multiple  access  and 
intentional  interference  (jamming). 

1.3  Research  Assumptions 

The  following  assumptions  were  used  when  modeling  the  communication  sys¬ 
tem  and  are  consistent  with  previous  TDCS  work  [7,  5,  2] 

1.  Communication  channel  noise  is  additive  white  Gaussian  noise. 

2.  Perfect  synchronization  between  each  transmitter /receiver  pair. 

3.  Perfect  code  synchronization  between  each  transmitter /receiver  pair. 

4.  There  are  no  multi-path  signals  present  in  the  network.  Transmitted  signals  only 
propagate  along  one  line-of-sight  path  to  the  receiver  of  interest. 

5.  There  is  no  Doppler  shift  on  the  received  signals.  The  receiver  of  interest  and  all 
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transmitters  are  stationary  relative  to  each  other. 


1.4  Research  Scope 

The  scope  of  this  research  was  limited  by  assumptions  contained  in  Section  1.3. 
Spectrally  encoded  TDCS  techniques  were  analyzed  and  simulations  conducted  for 
communication  performance  and  multiple  access  capability  using  M  =  2,4,8  and 
16  communication  symbols,  i.e.,  TD-BPSK,  TD-QPSK,  TD-8PSK,  and  TD-16PSK 
systems,  respectively.  The  first  research  goal  was  to  establish  the  probability  of 
symbol  error  (Pe)  and  probability  of  bit  error  (Pb)  versus  normalized  signal-to- 
noise  ratio  (Ej-JNq)  performance  for  the  TD-MPSK  systems  considered.  The  error 
was  determined  for  different  multiple  access  and  jamming  environments.  The  second 
research  goal  was  to  determine  Pb  performance  versus  the  number  of  multiple  access 
users  in  a  jamming  environment.  The  third  research  goal  involved  an  investigation 
for  optimizing  the  spectral  encoding  (phase  code  selection  and/or  assignment)  such 
that  multiple  access  interference  is  minimized. 

1.5  Research  Approach 

The  first  part  of  the  research  involved  establishing  communication  performance 
of  a  single  TD-MPSK  transmitter-receiver  pair.  Performance  for  this  part  of  the  re¬ 
search  was  characterized  using  probability  of  bit  error  (Pb)  versus  the  normalized 
signal-to- noise  ratio.  A  Matlab®  model  of  the  TD-MPSK  communication  system  was 
developed.  This  model  generated  a  random  data  modulated  signal,  added  AWGN 
to  the  signal  and  then  demodulated  the  signal.  After  demodulating  the  signal,  com¬ 
munication  symbol  estimates  were  compared  to  originally  modulated  symbols  to 
determine  if  there  was  an  error.  A  Monte  Carlo  analysis  was  then  conducted  using 
this  modulation-demodulation  process  which  was  repeated  at  different  normalized 
signal-to-noise  ratios.  The  Pb  for  different  normalized  signal-to-noise  ratios  was  then 
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observed.  The  observed  Pb  was  then  compared  with  that  obtained  from  expressions 
for  analytic  Pb  performance.  Different  scenarios  were  tested  for  multiple  access  in¬ 
terference,  jamming,  and  notching  of  frequency  components  with  interference.  The 
second  part  of  the  research  involved  determining  communication  performance  of  a 
single  transmitter-receiver  pair  using  a  fixed,  normalized  signal-to-noise  ratio  while 
varying  the  number  of  network  users.  Using  the  same  Matlab®  model  used  for  com¬ 
munication  performance  characterization,  random  symbols  were  again  generated,  the 
fundamental  signaling  waveform  was  modulated,  AWGN  was  added,  and  multiple 
access  signals  from  the  different  users  were  added.  The  demodulated  symbols  at  the 
receiver  were  then  compared  to  modulated  symbols  and  errors  determined.  A  Monte 
Carlo  analysis  was  again  conducted  using  the  modulation-demodulation  process  us¬ 
ing  a  different  numbers  of  additional  users.  The  Pb  for  the  different  number  of  users 
was  observed  and  once  again  compared  to  analytic  Pb  performance.  The  different 
scenarios  included  jamming,  notching,  and  a  combination  of  jamming  and  notching. 

1.6  Materials  and  Equipment 

The  proposed,  spectrally  encoded  TDCS  system  was  simulated  using  Matlab®, 
Version  6.1.  Mathworks,  Inc.,  developed  by  Matlab®  as  a  computer  language  for 
technical  computing [8].  The  simulations  were  run  on  a  1.0  GHz  Pentium  111  personal 
computer. 

1.7  Thesis  Organization 

Chapter  2  of  the  thesis  provides  background  information  on  the  conceptual 
TDCS  system,  phase  shift  keying  modulation,  CSK  modulation,  and  linear  phase 
shift  registers.  Chapter  3  provides  an  overview  of  how  the  communication  signals 
and  the  RF  environment  are  modeled.  It  also  provides  a  discussion  on  Monte  Carlo 
analysis  as  used  in  computer  simulation  of  the  proposed  TDCS  technique.  Chapter  4 
provides  a  comparison  of  theoretical  and  simulated  results  for  different  scenarios 
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considered  under  the  research.  Chapter  5  discusses  conclusions  that  can  be  drawn 
from  the  research  and  provides  recommendations  for  future  research.  Appendix  A 
contains  the  Matlab®  simulation  code  developed  and  used  for  the  communication 
system.  Appendix  B  contains  mathematical  derivations  for  some  of  the  key  results 
obtained  as  part  of  the  research  and  dispersed  throughout  the  document. 
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2.  Literature  Review 


2.1  Introduction 

This  chapter  introduces  MPSK  as  a  modulation  technique  obtained  via  spectral 
encoding,  with  specific  application  to  a  multiple  access  transform  domain  commu¬ 
nication  system  (MA-TDCS).  Section  2.2  describes  MPSK.  Section  2.3  describes  a 
conceptual  transform  domain  communication  system  (TDCS)  which  uses  spectral 
encoding  to  achieve  multiple  access  capability  and  interference  suppression.  Sec¬ 
tion  2.4  describes  code  division  multiple  access  (CDMA)  and  Section  2.5  summarizes 
the  chapter. 

2.2  Phase  Shift  Keying 

Phase  shift  keying  (PSK)  is  a  method  of  modulation  where  each  communica¬ 
tion  symbol  has  equal  energy  and  can  be  represented  as  the  sum  of  two  amplitude 
modulated  orthonormal  signals  (a  pair  of  unit  energy,  orthogonal  basis  functions). 
The  orthonormal  signals  are  represented  by  the  functions  and  t/bW-  The  sum 
of  the  two  amplitude  modulated  orthonormal  signals  making  up  the  ith  symbol  (one 
of  M )  is  given  by  (2.1)  [3].  The  symbol  phase  values  4>i  can  be  assigned  per  (2.2) 
which  ensures  equal  angular  spacing  between  phase  values. 

Si(t)  —  VEcos  +  VEsin  {4>i)  "02 (t)  0<f<T  (2.1) 

2m 

O;  u  i  =  1, M  (2.2) 

2.2.1  Binary  Phase  Shift  Keying.  For  Binary  phase  shift  keying 
(BPSK)  there  are  two  communication  symbols  with  each  representing  one  of  two 
bit  values.  Given  M  —  2  for  binary  signaling,  (2.1)  reduces  to  (2.3). 
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(2.3) 


Si(t)  =  \[E  cos  (  (pi)  "01  (t)  0  <  t  <  T  i  =  1,2 

4>i  —  ni  i  =  1,2  (2.4) 

For  optimal  maximum  likelihood  detection  and  estimation  of  BPSK  signals, 
only  one  correlator  is  needed  to  obtain  matched  filter  performance.  Assuming  the 
transmitted  bit  values  are  equally  probable,  and  signaling  is  being  conducted  over 
and  AWGN  channel,  the  single  correlator  can  be  “matched”  to  the  basis  function 
used  to  create  the  communication  symbols  and  achieve  optimal  performance  (min¬ 
imum  bit  error  rate).  In  this  case,  the  optimum  detection  threshold  is  zero  and 
the  correlator  output  “sign”  determines  the  symbol  estimate,  e.g.,  a  positive  output 
dictates  estimating  the  bit  as  a  1  and  a  negative  output  dictates  estimating  the  bit 
as  a  0.  Figure  2.1  shows  a  block  diagram  of  a  single  channel  BPSK  correlation  re¬ 
ceiver  [3]  where  the  correlation  reference  signal  is  the  ipi  (f)  basis  function.  Received 
signal  r(t)  is  the  sum  of  the  modulated  communication  signal,  environmental  inter¬ 
ference  (other  signals)  and  channel  noise.  Received  signal  r(t )  and  the  correlator 
reference  signal  are  multiplied  and  integrated  over  one  symbol  interval  Ts.  The 
sign  of  the  integration  result  forms  the  basis  for  making  a  decision  (estimating)  as 
to  which  of  the  two  symbols  are  present  in  r(t). 


Figure  2.1  Single  Channel  BPSK  Correlation  Receiver 

Given  there  are  only  two  symbols  in  binary  signaling,  and  assuming  equal 
energy  signaling,  the  average  symbol  energy  (Es)  equals  average  bit  energy  (Eb). 
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Once  the  normalized  signal-to-noise  ratio  [E^/Nq)  is  calculated,  probability  of  bit 
error  Pb  can  be  determined  using  (2.5)  [3]  where  Nq/2  is  the  two-sided  noise  power 
spectral  density. 


2.2.2  M-Ary  Phase  Shift  Keying.  For  M-Ary  Phase  Shift  Keying 
(MPSK),  there  are  M  communication  symbols  where  M  =  2k  and  k  >  1;  each 
MPSK  symbol  represents  k  bits.  Figure  2.2  shows  a  two  channel  correlation  receiver 
that  can  be  used  for  all  MPSK  systems  [3].  Received  signal  r(t)  is  split  and  passed 
to  two  correlation  paths.  In  each  correlation  path,  r(t)  is  multiplied  by  one  of 
the  two  basis  functions  originally  used  to  represent  the  communication  symbols. 
Following  integration  over  each  symbol  interval,  a  ratio  (Y / X)  is  formed  between 
the  two  correlator  channel  outputs.  Estimated  phase  value  (j)  is  established  by  taking 
the  arctangent  of  Y/X  and  compared  with  each  possible  phase  value  of  the  M 
communications  symbols.  The  symbol  phase  value  which  is  closest  to  the  estimated 
phase  value  is  estimated  as  being  received.  Under  equally  probable  conditions  (all 
symbols  transmitted  with  equal  apriori  probability)  and  signaling  over  an  AWGN 
channel,  this  two  channel  correlation  process  is  equivalent  to  maximum  likelihood 
detection  and  provides  minimum  probability  of  symbol  error  (Pe)  and  bit  error 
(Pb)  [9] 

For  MPSK  signaling,  average  symbol  energy  ( Es )  is  k  times  average  bit  energy 
(Eb)  where  k  is  the  number  of  bits  represented  by  each  symbol.  Equation  (2.6)  and 
(2.7)  provide  good  estimates  of  Pe  and  Pb  under  high  signal-to-noise  conditions  and 
are  based  on  Gray  code  bit-to-symbol  assignment  [3].  These  equations  are  primarily 
based  on  accounting  for  adjacent  symbol  errors  only,  i.e.,  they  only  take  into  account 
the  most  probable  (dominant)  symbol  error  conditions  which  are  adjacent  symbol 
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Figure  2.2  Two  Channel  Correlation  Receiver  for  MPSK  Signaling 


estimation  errors.  As  the  signal-to-noise  ratio  decreases,  adjacent  symbol  errors 
become  less  dominant  and  these  equations  become  less  accurate. 


Under  lower  signal-to-noise  conditions,  (2.8)  provides  a  more  accurate  means 
for  calculating  a  theoretical  Pb  estimate  [10]. 


-t  m  l  f  E  \ 

=  F  <2'8) 

m= 1  '  ' 

The  Dm  term  in  (2.8)  is  calculated  using  (2.9)  where  Q\x\  is  the  complementary 
error  function.  The  Wm  terms  in  (2.8)  are  given  in  Table  2.1  for  the  three  M  values 
shown. 
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Table  2.1  Mapping  Weights  for  Theoretical  Pg  Calculations  per  (2.8) 


M  =  4 

M  =  8 

M  =  16 

m 

wm 

nr 

wm 

m 

wm 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

3 

1 

3 

2 

3 

2 

4 

2 

4 

2 

5 

2 

5 

2.5 

6 

2 

6 

3 

7 

1 

7 

2.5 

8 

2 

9 

2.5 

10 

3 

11 

2.5 

12 

2 

13 

2 

14 

2 

15 

1 

Jo00/  (z-  {Q[z  tan((2m  —  1)^)]  -  Q[2;tan((2m  +  1)^)]}^ 

m  =  0,l,...,f  -  1 

f0-  f  (z+  \l2kNo)  {Qlz  tan((M  -  2m  -  1)^)]  -  Q[ztan((M  -  2m  +  1)^ )}}  dz 
m=  f  + 

(2.9) 


/(<)  =  e  2* 

It  should  be  noted  that  (2.9)  in  and  of  itself  does  not  provide  all  the  required 
values  of  Dm  for  all  m.  However,  the  relationships  provided  in  (2.10)  and  (2.11)  can 
be  used  to  permit  complete  generation  of  required  Dm  values. 
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Figure  2.3  Conceptual  TDCS  Transmitter 


DM-m  =  An  (2.10) 

M—l 

=  1  (2-11) 

m= 0 

2.3  Transform  Domain  Communication  System  (TDCS) 

The  purpose  of  a  TDCS  system  is  to  communicate  in  a  radio  frequency  (RF) 
environment  containing  extensive  interference.  Normally,  communication  systems 
attempt  to  mitigate  interference  effects  by  employing  sophisticated  filtering  tech¬ 
niques  in  the  time  and/or  frequency  domains.  A  TDCS  effectively  samples  the  RF 
environment  and  generates  “smart”  waveforms  at  the  transmitter;  “smart”  here  sim¬ 
ply  implies  the  transmitted  waveforms  contain  no  (minimal)  energy  in  spectral  areas 
containing  interference  [5].  Figure  2.3  shows  a  block  diagram  of  a  conceptual  TDCS 
transmitter. 
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2.3.1  TDCS  Transmitter.  The  TDCS  transmitter  of  Fig.  2.3  uses 
several  processing  steps  for  signal  generation  prior  to  transmission.  First,  the  RF 
environment  is  sampled  and  a  spectral  estimate  of  the  energy  obtained.  Second,  the 
spectral  estimate  is  analyzed  to  determine  which  region(s)  of  the  RF  environment 
contain  interference  (interference  here  includes  all  sources  of  radiation  not  generated 
by  the  TDCS  itself).  Third,  a  fundamental  signaling  waveform  is  generated  con¬ 
taining  frequency  components  in  the  region (s)  containing  interference.  Fourth,  the 
fundamental  waveform  is  stored  for  subsequent  replication  and  modulation.  Finally, 
the  fundamental  waveform  is  data  modulated  to  create  the  final  signal  for  transmis¬ 
sion  [5].  The  following  paragraphs  describe  the  functionality  of  each  block  in  Fig.  2.3 
in  greater  detail. 

2. 3. 1.1  Spectral  Estimation.  The  first  step  in  TDCS  processing 
involves  estimating  /  locating  where  in  frequency  spectral  interference  is  occurring. 
In  fourier-based  spectral  estimation,  the  RF  environment  is  sampled  at  uniform 
intervals  and  a  discrete  Fourier  transform  (DFT)  performed.  The  DFT  is  a  com¬ 
putational  tool  for  performing  frequency  analysis  of  a  sampled  time  domain  signal. 
The  DFT  effectively  describes  the  spectral  content  of  the  RF  environment  using  a 
discrete  number  of  equally  spaced  frequency  components.  The  number  of  frequency 
components  equals  the  number  of  samples  taken.  Each  frequency  component  is  de¬ 
scribed  by  frequency,  amplitude,  and  phase  parameters  [11].  The  TDCS  processing 
extracts  frequency  and  amplitude  information  from  the  DFT  and  determines  which 
spectral  region(s)  (actual  frequency  components)  is  clear  of  interference. 

2.3. 1.2  Spectral  Thresholding.  Given  the  spectral  magnitude 
response  provided  by  the  DFT  processing,  the  TDCS  establishes  a  threshold  such 
that  all  frequency  components  (samples)  exceeding  the  threshold  (those  containing 
interference)  are  assigned  a  value  of  zero.  Frequency  components  whose  magnitudes 
are  below  the  threshold  are  deemed  interference  free  and  assigned  a  value  of  one  [12], 
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2.3. 1.3  Spectral  Phase  Coding.  Given  the  TDCS  thresholding 
process  has  established  which  frequency  components  are  available  for  fundamental 
waveform  generation  (dictated  by  the  DFT  sample  values  that  have  been  set  to  one), 
selection  of  phase  values  to  be  assigned  to  those  components  provides  the  degree  of 
freedom  required  to  enable  interference  avoiding,  multiple  access  capability;  the  se¬ 
lection  of  phase  values  and  their  subsequent  assignment  to  frequency  components  is 
the  form  of  spectral  encoding  employed  here.  The  Random  Phase  Generator  gen¬ 
erates  a  phase  for  each  of  these  frequency  components.  In  previous  TDCS  work,  a 
Pseudo-Random  Phase  Generator  has  been  employed  to  generate  the  required  sets 
of  phase  values  [6].  In  this  case,  linear  feedback  shift  registers  were  used  to  gener¬ 
ate  periodic  pseudo-random  sequences  [4]  which  were  subsequently  used  in  a  “phase 
mapping”  process  which  pseudorandomly  assigned  a  unique  sequence  of  phase  values 
to  each  user  in  the  network.  If  both  the  transmitter  and  receiver  use  the  same  initial 
register  contents  for  their  linear  shift  register,  and  the  receiver  can  autonomously 
synchronize  to  the  transmitter,  reliable  detection  and  estimation  takes  place  [6]. 

2.3. 1.4  Magnitude  Scaling  for  Equal  Energy  Signaling.  After 
spectral  estimation,  thresholding  and  spectral  encoding,  the  frequency  components 
are  scaled  to  ensure  that  equal  energy  symbols  are  transmitted  independent  of  the 
number  frequency  components  being  used.  For  example,  if  only  a  third  of  the  fre¬ 
quency  components  remain  after  thresholding  due  to  interference,  the  amount  of 
required  power  contained  in  each  remaining  component  is  tripled  relative  to  the  case 
when  no  interference  is  present. 

2.3. 1.5  Inverse  Discrete  Fourier  Transform.  The  fundamen¬ 
tal  signaling  waveform  is  obtained  by  taking  the  inverse  discrete  Fourier  transform 
(IDFT)  of  the  scaled,  spectrally  encoded  frequency  components.  The  energy  in 
the  resultant  time  domain  waveform  is  contained  entirely  in  frequency  components 
(spectral  locations)  where  no  interference  is  present  [12],  The  final  set  of  M  commu- 
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nication  symbols  is  then  created  from  the  fundamental  signaling  waveform.  Perhaps 
the  easiest  form  of  data  modulation  that  can  be  employed  is  binary  antipodal  sig¬ 
naling.  For  binary  antipodal  signaling,  one  symbol  is  the  fundamental  waveform 
itself  and  the  other  symbol  is  the  negative  of  the  fundamental  waveform.  If  only 
one  frequency  component  remains  after  thresholding,  the  fundamental  waveform  is 
a  simple  sinusoid  and  antipodal  signaling  produces  results  identical  to  BPSK  [3]; 
multiplying  a  simple  sinusoid  by  negative  one  is  equivalent  to  shifting  its  phase  by  n. 
Using  the  TDCS  process  described  above,  it  can  be  shown  that  antipodal  signaling 
can  be  achieved  by  either  1)  negating  the  fundamental  waveform  after  generation,  or, 
2)  by  shifting  the  spectrally  encoded  phases  of  all  frequency  components  by  7r.  An 
alternate  form  of  data  modulation  that  has  been  explored  extensively  in  previous 
TDCS  applications  is  CSK  [5,  2,  6].  In  CSK  symbols  are  generated  by  cyclically 
shifting  (in  the  time  domain)  the  fundamental  signaling  waveform  in  proportional 
amounts  determined  by  the  number  of  desired  communication  symbols.  For  exam¬ 
ple,  in  binary  CSK  the  first  symbol  would  be  the  fundamental  signaling  waveform 
itself  and  the  second  symbol  would  be  generated  by  extracting  the  second  half  of  the 
fundamental  waveform  and  concatenating  it  onto  the  first  half  of  the  waveform. 

2.3. 1.6  Waveform  Storage  in  Memory.  After  the  fundamental 
signaling  waveform  is  generated,  it  is  stored  in  memory  and  made  available  to  the 
modulator  for  subsequent  data  modulation.  As  changes  in  the  RF  environment  are 
sensed,  or  at  scheduled  intervals,  the  fundamental  waveform  generation  process  is 
repeated  and  memory  updated  [12]. 

2.3. 1.7  Data  Modulation.  Data  (information)  is  modulated  onto 
the  fundamental  waveform  prior  to  transmission.  Data  modulation  is  the  process 
whereby  baseband  information  is  impressed  onto  a  bandpass  signal  [1],  In  the  TDCS 
case,  the  baseband  information  (data  bits)  is  received  by  the  modulator  and  the 
modulator  decides/determines  which  of  the  M  communication  symbols  to  output 
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in  accordance  with  a  predefined  bit-to-symbol  assignment  process.  For  the  binary 
antipodal  signaling  described  earlier,  the  modulator  would  simply  take  in  one  bit  at 
a  time  and  output  either  Si(f)  (for  a  binary  1)  or  s2(t)  =  —  Si(t)  (for  a  binary  0). 

2.3. 1.8  Signal  Transmission.  Following  data  modulation,  the 
continuous  digitally  encoded  waveform  is  ready  for  transmission.  The  transmitter 
may  do  several  things  to  the  waveform  prior  to  actual  transmission,  including,  1)  ap¬ 
ply  amplification/gain,  2)  apply  shaping/filtering,  3)  apply  frequency  translation  for 
radiation  efficiency,  4)  apply  frequency  translation  to  ensure  operation  in  assigned 
spectral  region,  etc. 

2.3.2  TDCS  Receiver  Processor.  Functionally,  the  TDCS  receiver 
contains  the  same  fundamental  waveform  generator  as  the  transmitter  and  employs 
a  conventional  M- channel  correlation  receiver  for  symbol  detection  and  estimation. 
Overall  system  performance  is  primarily  driven  by  how  “similar”  the  fundamental 
waveform  generated  in  the  receiver  “looks  like”  the  one  generated  at  the  transmitter. 
Under  ideal  conditions,  the  receiver  and  transmitter  “see”  the  same  interference  and 
generate  identical  fundamental  waveforms;  the  transmitter  puts  no  energy  in  the 
frequency  components  that  were  notched  out  and  the  receiver  extracts  no  energy 
(interference)  from  these  same  components.  In  this  case,  the  only  remaining  error 
is  attributable  to  the  AWGN  channel  itself.  Assuming  each  of  the  TDCS  receiver 
channels  is  “matched”  to  one  of  M  possible  communication  symbols,  maximum 
likelihood  detection  and  estimation  is  accomplished  by  simply  estimating  the  received 
symbol  based  on  the  maximum  correlator  output  [3]. 

2.4  Code  Division  Multiple  Access 

In  code  division  multiple  access  (CDMA),  multiple  users  share  communication 
resources  (e.g.,  spectrum)  without  significantly  affecting  each  other’s  performance. 
In  general,  multiple  users  in  a  CDMA  systems  do  not  need  to  be  synchronized  to  ef- 
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fectively  operate.  Coding  of  the  multiple  users  is  usually  done  in  the  time  domain  [4]. 
As  demonstrated  in  previous  TDCS  research,  it  is  possible  to  apply  phase  coding  in 
the  frequency  domain  (spectrally  encode)  and  achieve  multiple  access  performance; 
by  selecting  a  set  of  codes  with  desirable  properties  and  uniquely  assigning  them  to 
users  CDMA  capability  can  be  achieved  with  TDCS  techniques  [6]. 

Frequency  Division  Multiple  Access  (FDMA)  and  Time  Division  Multiple  Ac¬ 
cess  (TDMA)  are  limited  in  their  multiple  access  capacity  by  bandwidth.  The  mech¬ 
anism^)  involved  in  implementing  CDMA  is  different  from  these  and  multiple  access 
capability  is  limited  by  interference  [13].  The  capacity  is  limited  only  by  the  level  of 
interference  the  communication  system  can  tolerate. 

2.5  Summary 

The  conceptual  TDCS  employing  spectral  encoding,  as  introduced  in  this  chap¬ 
ter,  provides  one  method  for  mitigating  interference  affects  while  achieving  multiple 
access  communications.  By  creating  “smart”  waveforms  at  the  transmitter  that  ef¬ 
fectively  “avoid”  spectral  regions  containing  interference,  the  receiver’s  task  is  greatly 
simplified  in  that  sophisticated  filtering  techniques  are  not  required  to  remove  inter¬ 
ference.  The  use  of  MPSK  with  TDCS  process  promises  to  provide  more  efficient 
use  of  bandwidth  relative  to  other  modulation  techniques  (antipodal  signaling  and 
CSK  modulation)  considered  previously  for  TDCS  applications. 
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3.  Modeling  and  Simulation  Techniques 

3.1  Introduction 

This  chapter  provides  an  overview  on  how  TDCS  symbols  interact  with  each 
other  and  with  interference.  It  also  describes  how  the  TDCS  can  be  used  in  a  MPSK 
system.  Simulation  assumptions  and  a  description  of  the  Monte  Carlo  simulation 
process  are  also  discussed. 

3.2  Transform  Domain  Symbol  Definition 

Symbol  development  for  the  proposed  TD-MPSK  technique  begins  in  the  fre¬ 
quency  domain  where  phase  modulation,  both  data  and  coded  multiple  access,  are 
applied  directly  to  sinusoidal  components;  this  is  the  spectral  encoding  referenced 
throughout  this  document.  Although  the  spectral  encoding  process  used  here  is  sim¬ 
ilar  to  orthogonal  frequency  division  multiplexing  (OFDM)  [14],  i.e.,  in  the  OFDM 
technique  data  phase  modulation  is  applied  in  the  frequency  domain,  the  TD-MPSK 
spectral  encoding  process  differs  in  that  data  phase  modulation  applied  here  is  con¬ 
stant  across  all  spectral  components  for  a  given  user.  The  kth  TD-MPSK  communi¬ 
cation  symbol  for  user  v  is  generated  from  (3.1)  where  the  frequency,  time  and  phase 
parameters  are  defined  in  Table  3.1.  The  time  domain  representation  of  a  TD-MPSK 
symbols  is  generated  by  taking  the  inverse  Fourier  transform  of  (3.1),  as  shown  in 
(3.2),  and  selecting  a  portion  of  this  result  to  represent  the  actual  communication 
symbols  as  shown  in  (3.3). 
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Table  3.1  TD-MSPK  Symbol  Parameters 


t 

Time 

p 

Number  of  Sinusoidal  Components 

Tsb 

Symbol  Duration 

fsb 

Spectral  Sample  Spacing 

W 

Data  Phase  Modulation 

~W 

Multiple  Access  Phase  Coding 

p 

(t)  =  2  cos  (271  fsbpt  +  ^  to  <t  <to  +  Tsb  (3.3) 

p=i 

The  time  domain  symbol  sampling  rate  is  chosen  as  an  integer  multiple  of 
fsb  where  sampling  starts  at  t  —  tQ\  sampling  can  start  at  t0  —  0  without  loss 
of  generality.  The  important  information  in  the  symbol  is  amplitude,  frequency 
and  phase  for  each  sinusoid  in  the  symbol.  To  ensure  all  the  information  from  the 
continuous  waveform  is  in  the  sampled  waveform,  the  sampling  needs  to  be  greater 
than  twice  the  frequency  of  the  highest  frequency  sinusoid.  The  equation  for  time 
sampled  TD-MSPK  symbols  is  given  in  (3.4). 
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The  discrete  fourier  transform  (frequency  index  =  m)  of  the  symbol  is  shown 
in  equation  (3.5) 
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Since  each  TD-MPSK  symbol  is  sampled  at  a  frequency  which  is  more  than 
twice  that  of  the  highest  frequency  sinusoid,  IV  —  p  is  always  greater  than  p  and 
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there  is  no  aliasing.  Analysis  of  (3.5)  provides  the  key  as  to  how  time  domain  phase 
modulation  is  achieved  through  spectral  encoding  in  the  frequency  domain. 

The  energy  in  a  symbol  is  determined  from  its  autocorrelation.  The  energy  in 
any  TD-MPSK  symbol  is  shown  in  (3.6).  This  equation  is  derived  from  Section  B.3. 
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3.3  Symbol  Cross-Correlation  For  Synchronous  Users 

Modifying  (3.4),  a  symbol  from  another  user  is  shown  in  equation  (3.7). 
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The  cross-correlation  between  symbols  of  synchronous  users  is  shown  in  (3.8). 
This  equation  is  derived  from  Section  B.l. 
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3.3.1  Mean  and  Variance  of  Symbol  Cross-Correlation.  If  the 

data  modulation  is  completely  random,  the  mean  value  of  the  cross-correlation  be¬ 
tween  symbols  of  synchronous  users  is  zero.  In  this  case,  the  variance  of  the  cross¬ 
correlation  between  symbols  of  synchronous  users  for  the  TD-BPSK  version  is  given 
by  (3.9).  This  equation  is  derived  from  Section  B.4.  For  TD-MPSK  (M  =  2k  where 
k  >  1),  the  variance  of  the  cross  correlation  is  given  by  (3.10).  This  equation  is 
derived  from  Section  B.5. 
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If  the  data  phase  modulation  and  multiple  access  phase  coding  are  completely 
random,  the  mean  value  of  the  cross-correlation  between  symbols  of  synchronous 
users  is  zero.  The  variance  of  the  cross-correlation  between  symbols  of  synchronous 
users  for  this  case  is  given  by  (3.11).  This  equation  is  derived  from  Section  B.6. 


3.3.2  Mean  and  Variance  of  Symbol  Cross-Correlation:  Orthogo¬ 
nal  Code  Selection.  The  cross-correlation  between  communication  symbols  of 
different  synchronous  users  can  be  made  zero  through  proper  multiple  access  phase 
code  selection  (the  orthogonal  code  case).  Section  B.7  provides  a  derivation  of  how 
the  phase  codes  are  selected.  Orthogonal  symbols  are  defined  using  (3.12).  P  is  the 
number  of  sinusoids  that  make  up  a  symbol.  N  is  the  number  of  samples  in  a  symbol 
where  N  >  2 P.  v  is  the  number  of  the  user  where  v  can  equal  0  to  G  —  1  where 
G  is  defined  in  (3.12).  A  symbol  for  another  user  can  be  defined  with  (3.13)  where 
w  is  the  number  of  users  which  can  also  equal  0  to  G  —  1.  As  long  as  v  ^  w  the 
correlation  between  the  symbols  of  the  two  different  users  is  0. 
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3.4  Symbol  Cross-Correlation  and  Bandlimited  Noise 

Equation  (3.14)  provides  a  definition  for  narrowband  interference  (Inb)  con¬ 
sidered  for  this  work.  The  interference  is  modeled  similar  to  a  symbol  in  the  sense 
that  it  is  the  sum  of  sinusoids,  each  with  a  specific  amplitude.  However,  the  phases 
of  the  sinusoidal  components  making  up  the  narrow  band  interference  change  each 
time  the  interference  is  generated. 


Inb  (n)  =  —  Bp  cos  (2vr ^  (3.14) 

p= i 

If  the  amplitude  of  each  sinusoidal  component,  Bpi  is  known,  the  noise  can  be 
treated  as  a  synchronous  symbol  with  random  spreading  phases.  This  gives  the  cross¬ 
correlation  mean  and  variance  results  obtained  when  cross-correlating  synchronous 
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symbols  with  random  multiple  access  phase  coding.  The  final  cross-correlation  result 
is  shown  in  (3.15). 


(3.15) 


3.5  Symbol  Cross-Correlation  for  Asynchronous  Users 

For  the  asynchronous  case  three  symbols  need  to  defined.  The  first  symbol 
comes  from  user  v  and  is  defined  by  (3.16).  The  second  and  third  symbols  come 
from  user  w  and  are  defined  by  (3.17)  and  (3.18)  respectively.  The  first  symbol  is 
correlated  with  the  concatenation  of  the  last  part  of  the  second  symbol  and  the  first 
part  of  the  third  symbol.  The  end  part  of  the  second  symbol  has  a  length  of  a  symbol 
interval  minus  a  delay  defined  as  r.  The  beginning  part  of  the  third  symbol  has  a 
length  of  r. 
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Using  the  definitions  of  (3.16),  (3.17)  and  (3.18),  the  cross-correlation  between 
Symbol  i  of  User  v  and  the  symbol  transition  boundary  of  Symbol  j  and  Symbol  k 
for  User  w,  can  be  expressed  as 
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For  the  case  of  the  sampled  symbol  where  multiple  access  phase  coding  is  ran¬ 
dom  and  uniformally  distributed  between  [0,  27t],  the  variance  of  the  cross-correlation 
can  be  estimated  with  (3.20). 
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When  the  waveforms  are  sampled  at  a  rate  of  N  samples  per  symbol  where 
P  <  2 N,  the  correlation  can  be  estimated  by  substituting  N  for  fsb.  The  estimate 
of  the  correlation  of  the  sampled  waveforms  is  given  by  3.21. 
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3.6  Cross-Correlation  Between  Two  Symbols  of  Single  Synchronous  User 

As  derived  in  Appendix  B,  the  cross-correlation  between  two  symbols  of  a 
single  synchronous  user  is  shown  in  (3.22).  Using  the  symbol  energy  expression 
of  (3.6),  (3.22)  can  be  rewritten  as  shown  in  (3.23).  As  indicated  in  (3.23),  the 
symbol  cross-correlation  is  simply  the  symbol  energy  multiplied  by  the  cosine  of  the 
difference  between  data  phase  modulation  values  for  the  two  symbol  intervals  under 
consideration. 


fl.(r>  =  4l cos  («W -»<”>)  (3.22) 

p= 1 

RA-u)  =  Es  cos  (e\v)  -  (3.23) 

The  cross-correlation  results  of  (3.23)  lead  directly  into  the  MPSK  data  modu¬ 
lation  technique  proposed  for  the  TDCS.  The  first  step  in  developing  a  MPSK  system 
is  determining  two  orthogonal  symbols.  Equation  (3.24)  shows  the  equation  for  the 
first  orthogonal  symbol.  Equation  (3.25)  shows  the  equation  for  the  second  orthogo¬ 
nal  symbol.  These  two  equations  are  the  reference  symbols  used  in  the  demodulator 
of  the  MPSK  system. 


^{v)  =  ^J2ap  cos  {27Tpjv + 4w))  (3-24) 

p= i 

^2V)  =  ^  AP  sin  (2np^  +  4W))  (3-25) 

p= 1 

The  correlation  with  the  hrst  symbol  is  the  symbol  energy  multiplied  by  the 
cosine  of  the  angle  modulation.  The  correlation  with  the  second  symbol  is  the  symbol 
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energy  multiplied  by  the  sine  of  the  angle  modulation.  The  arctangent  of  the  second 
correlation  divided  by  the  first  correlation  is  the  phase  modulation. 


3.7  Interference  Model 

If  the  normalized  signal-to-noise  ratio  is  known,  the  probability  of  symbol 
error  (Pe)  and  probability  of  bit  error  (Pg)  can  be  calculated  in  the  presence  of 
interference.  When  multiple  access  interferers  are  present,  the  normalized  signal- 
to-noise  ratio  can  be  modified  to  account  for  the  additional  interference  effects.  In 
this  case,  the  Eb/N0  term  in  analytic  Pe  and  Pg  expressions  of  section  2.2.1  and 
section  2.2.2  are  commonly  replaced  with  E^/^Nq  +  Ni)  to  incorporate  interference 
effects;  this  Eb/(N0  +  W)  can  be  rewritten  as  shown  in  (3.26)  [6]. 


Eb  (ea  1  (Eb  y1 

No  +  Ni  \N0 )  \Ni) 


(3.26) 


For  the  proposed  TD-MPSK  technique,  the  Eb/Ni  term  of  (3.26)  can  be  es¬ 
timated  using  (3.27)  where  R  is  the  number  of  bits  per  communication  symbol. 
£e4>  |  ^s(«)swj  |  is  calculated  using  (3.11).  Es  is  calculated  using  (3.6). 


Eb 

Ni 


(3.27) 


If  the  multiple  access  phase  coding  (<pp)  is  known  for  both  users  the  Eb/Ni 
term  of  (3.26)can  be  estimated  using  (3.28).  £q  f  1  is  calculated  using 

(3.9)  for  the  TDCS-BPSK  case  and  (3.10)  for  the  TD-MPSK  case.  Es  is  once  again 
calculated  using  (3.6). 


Eb 

Ni 


Es 


°r.(«)  „(«>) 


(3.28) 
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3.8  Monte  Carlo  Verification 


To  verify  correctness  of  the  analytical  results  obtained  thus  far,  Monte  Carlo 
simulations  were  conducted  with  Matlab®  generating  the  AWGN  channel  response. 
The  AWGN  was  simply  summed  with  the  transmitted  signal  upon  reception  by  the 
receiver.  The  total  received  signal  was  demodulated  and  the  estimated  communica¬ 
tion  symbol  compared  with  the  modulated  symbol.  Multiple  Bernoulli  trials  were 
conducted  and  probability  of  error  determined  [7]. 

3.9  Modeling  and  Simulation  Assumptions 

Several  limitations  are  placed  on  communication  system  simulations.  First, 
spectral  estimation  was  deemed  outside  the  scope  of  this  research.  Therefore,  a 
predetermined  interference  envelope  was  established  such  that  the  spectral  notch 
removes  most  of  the  jamming  power.  Second,  phase  mapping  development  and 
characterization  was  also  deemed  outside  the  scope  of  this  thesis.  Extensive  research 
has  been  conducted  on  the  phase  mapping  process  used  with  the  conceptual  TDCS 
technique  [6];  these  previous  results  are  believed  to  be  directly  extendable  into  the 
proposed  TD-MPSK  process.  The  random  multiple  access  phase  codes  used  here 
are  generated  based  on  Matlab®’s  “RAND"  function.  Third,  transmitter /receiver 
pairs  are  synchronous.  Finally,  data  phase  modulation  values  are  equally  likely  and 
uniformly  distributed  for  all  M  values  considered. 

3.10  Summary 

This  chapter  provided  the  definition  and  development  of  proposed  a  TDCS 
technique  using  spectral  encoding  for  both  data  phase  and  multiple  access  phase 
modulations.  This  chapter  described  how  phase  values  were  selected  and  MPSK 
was  achieved  in  a  multiple  access  environment;  the  resultant  technique  is  termed 
TD-MPSK.  Communication  symbol  cross-correlation  was  characterized  in  both  syn¬ 
chronous  and  asynchronous  networks.  The  correlation  of  communication  symbols 
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with  interference  was  also  discussed.  Next  the  chapter  discussed  how  MPSK  sig¬ 
naling  can  be  implemented  using  conventional,  multichannel  correlation  techniques. 
Finally,  the  assumptions  used  for  modeling  and  simulation  were  presented. 
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4.  Modeling  and  Simulation:  Random  Code  Selection 

4.1  Introduction 

This  chapter  provides  an  overview  of  the  modeling  parameters  and  the  sce¬ 
narios  used  for  testing  the  proposed  TD-MPSK  system.  The  parameters  described 
include  the  communication  symbol  spectral  response,  the  jammer  spectral  response 
and  the  notched  symbol  spectral  response.  The  modeled  scenarios  include:  1)  no 
interference  present  (AWGN  only),  2)  multiple  access  interference,  3)  jamming  inter¬ 
ference  and  4)  multiple  access  and  jamming  interference  without  spectral  notching 
and  5)  multiple  access  and  jamming  interference  with  spectral  notching.  The  sce¬ 
narios  considered  help  characterize  the  TD-MPSK  system’s  ability  to  communicate 
in  the  presence  of  noise  and  interference.  All  scenarios  were  tested  with  and  without 
notching. 

4.2  Modeling  Parameters 

4.2.1  Spectral  Response  of  Communication  Symbol.  For  modeling 
and  simulation  purposes,  each  TD-MPSK  communication  symbol  is  made  up  of 
P  =  31  equal  amplitude  sinusoids.  The  sinusoid  frequencies  range  from  the  symbol 
rate  (fsb)  to  31  times  the  symbol  rate.  The  communication  symbols  are  sampled 
at  a  rate  of  64  samples  per  symbol.  By  sampling  the  symbols  at  this  rate,  there 
is  no  aliasing.  Equation  (4.1)  is  a  discrete  form  of  the  analytic  representation  of 
communication  symbols.  The  magnitude  of  the  discrete  Fourier  transform  of  each 
frequency  component  of  the  sampled  signal  equals  one  except  for  the  DC  and  Nyquist 
components  which  equal  zero.  Each  communication  symbol  is  made  up  of  64  samples. 
Since  the  number  of  samples  is  a  power  of  two,  processing  is  potentially  faster. 
Figure  4.1  is  a  stem  plot  showing  the  spectral  response  of  communication  symbols 
obtained  by  taking  the  magnitude  of  the  DFT  of  a  communication  symbol.  Figure  4.2 
is  a  corresponding  time  domain  representation  which  was  obtained  using  random 
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spectral  encoding,  i.e.,  a  random  data  phase  modulation  and  a  random  set  of  multiple 
access  phase  coding. 


-J>os  +  W  +  ek)) 


(4.1) 


Figure  4.1  Spectral  Response  (Magnitude  of  DFT)  of  Modeled  Communication 
Symbols 

4.2.2  Spectral  Response  of  Jamming.  Narrow  band  jamming  interfer¬ 
ence  can  be  represented  using  P  =  31  sinusoids  with  amplitudes  and  random  phase 
values  assigned  to  provide  desired  spectral  characteristics.  Figure  4.3  is  a  stem  plot 
of  the  magnitude  of  the  spectrum.  Equation  (4.2)  analytically  describes  the  jamming 
considered  using  the  values  of  Ap  listed  in  Table  4.1.  The  phase  values  assigned  to 
spectral  components  of  the  jammer  are  randomly  changed  for  every  iteration  of  the 
simulation. 


J(n) 


1 

32 


p 


p=31 


n 

64 


(4.2) 
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Figure  4.2  Time  Domain  Response  of  Modeled  Communication  Symbol 


Figure  4.3  Spectral  Response  of  Narrow  Band  Jammer  (P  =  31  Sinusoids) 

4.2.3  Spectrum  of  Notched  User.  As  stated  previously,  spectral  es¬ 
timation  was  not  within  the  scope  of  this  research.  However,  when  modeling  the 
spectrally  notched  cases,  frequency  components  of  the  communication  symbols  in 
regions  containing  the  highest  jamming  were  selectively  notched  out  (manually  set 
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Table  4.1  Amplitude  Coefficients  for  Jammer  of  Fig.  4.3 


Amplitude 

Value 

Amplitude 

Value 

Amplitude 

Value 

Ao 

0.0000 

A\ 

0.3034 

A2 

0.4142 

A3 

0.6362 

A4 

1.0888 

A5 

2.0215 

Aq 

3.4755 

a7 

3.8799 

As 

3.6811 

Ag 

2.9683 

A10 

1.9194 

An 

1.2279 

A±2 

0.8696 

A13 

0.6794 

A\4 

0.5667 

A15 

0.4928 

Ai6 

0.4406 

A17 

0.4017 

Ais 

0.3716 

Aig 

0.3477 

A20 

0.3284 

A-21 

0.3125 

A22 

0.2994 

A23 

0.2884 

A-24 

0.2793 

A25 

0.2718 

A26 

0.2656 

A27 

0.2605 

A28 

0.2566 

A29 

0.2535 

A30 

0.2514 

A31 

0.2502 

A32 

0.0000 

to  zero  after  applying  a  manual  threshold  level).  For  narrow  band  simulation,  8  of  31 
(~  25%)  sinusoidal  components  were  notched  out.  The  only  criteria  used  for  notch¬ 
ing  out  a  frequency  component  was  that  the  jamming  amplitude  for  that  particular 
component  exceeded  one.  Notching  out  these  frequency  components  significantly 
reduces  the  effects  of  the  jamming  interference  (12  dB  reduction  in  jamming  power). 
After  spectral  notching,  the  amplitude  of  remaining  frequency  components  is  scaled 
such  that  the  energy  of  notched  symbols  equals  that  of  the  original  symbols.  For 
this  25%  narrow  band  case,  the  amplitude  of  remaining  frequency  components  was 
raised  to  approximately  1.16.  Figure  4.4  is  a  stem  plot  of  the  notched  symbol. 

4.2.4  Number  of  Users.  One  variable  used  throughout  this  document 
is  Njj,  the  number  of  users  in  a  network.  For  these  simulations,  the  first  user  in 
the  network  is  defined  as  the  receiver  in  the  transmitter/receiver  pair  whose  com¬ 
munication  performance  is  being  analyzed.  The  second  user  in  the  network  is  the 
transmitter  in  the  transmitter /receiver  pair  being  analyzed.  Users  3  through  32  are 
additional  transmitters  in  the  network,  i.e.,  the  multiple  access  interferers.  When 
Njj  =  2,  there  is  no  multiple  access  interference.  When  Njj  =  3  the  number  of  trans- 
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Frequency  index  (k) 

Figure  4.4  Spectral  Response  (Magnitude  of  DFT)  of  Modeled  Communication 
Symbol  with  Narrow  Band  Notched  Spectrum  (P  =  31  Sinusoids) 

mitters  causing  multiple  access  interference  is  equal  to  one.  In  general,  the  number 
of  multiple  access  interferers  is  (iVj j  —  2). 

4.3  No  Interference  Present 

Figure  4.5  and  Fig.  4.6  show  probability  of  symbol  error  (Pe)  and  probability 
of  bit  error  (Pb)  versus  normalized  signal-to-noise  ratio  respectively.  The  Monte 
Carlo  results  matched  the  analytic  results  for  coherent  detection  of  TD-MPSK.  One 
thing  that  should  be  noted  is  that  Pb  and  Pe  are  the  same  for  TD-BPSK;  one  bit 
is  represented  by  one  symbol  in  TD-BPSK.  It  is  also  interesting  to  note  that  the 
Pb  versus  the  normalized  signal  to  noise  ratio  ( E^/Nq )  is  the  same  for  TD-BPSK 
and  TD-QPSK  [3].  Another  interesting  feature  of  TD-MPSK  is  that  as  the  level  of 
modulation  (M)  increases,  the  probability  of  error  increases  for  fixed  E^/Nq.  This 
means  that  there  is  a  trade-off  between  probability  of  error  and  bandwidth  efficiency. 


4-5 


E/rydE) 

Figure  4.5  No  Interference  Present:  Pe  vs  E^/Nq  for  TD-MPSK  Signaling  with 
P  =  31  Sinusoids 


Eb/N(j(dB) 


Figure  4.6  No  Interference  Present:  Pb  vs  E^/Nq  for  TD-MPSK  Signaling  with 
P  =  31  Sinusoids 

4.4  Spectral  Notching  Effects 

When  arbitrary  frequency  components  are  spectrally  notched  and  there  is  no 
interference  present,  the  Pe  and  Pb  performance  of  the  TD-MPSK  system  are  unaf- 
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fected.  Although  this  is  an  obvious  result,  given  that  scaling  is  applied  after  spectral 
notching,  simulations  were  conducted  to  ensure  proper  model  operation.  Figure  4.7 
and  Fig  4.8  show  PB  and  Pb,  respectively,  versus  normalized  signal-to-noise  ratio. 


E/rydS) 


Figure  4.7  Spectral  Notching:  Pe  vs  Eb/N0  for  TD-MPSK  Signaling  with  P  =  31 
Sinusoids 


WdB> 


Figure  4.8  Spectral  Notching:  PB  vs  Eb/N0  for  TD-MPSK  Signaling  with  P  =  31 
Sinusoids 
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4.5  Jamming  Present  With  No  Spectral  Notching 

For  the  following  simulations,  the  signal-to-interference  (average  symbol  energy- 
to-average  jamming  power)  ratio  ( J/S )  is  maintained  at  3.14  dB.  Figure  4.9  and 
Fig.  4.10  show  Pe  and  Pg  versus  normalized  signal-to-noise  ratio,  respectively,  for 
the  case  where  the  jammer  is  present.  Both  figures  show  how  error  performance  de¬ 
grades  when  interference  is  present.  Because  the  normalized  signal-to-noise  ratio  is 
maintained,  interference  effects  the  higher  levels  of  modulation  more  than  the  lower 
levels  of  modulation.  This  can  be  seen  by  comparing  Pe  for  a  TD-BPSK  system 
with  that  of  the  TD-16PSK  system. 


Figure  4.9  Jamming  Present  with  No  Spectral  Notching :  Pe  vs  E^/Nq  for  TD- 
MPSK  Signaling  with  J/S  =  3.14  dB  and  P  =  31  Sinusoids 


4.6  Jamming  Present  With  Spectral  Notching 

When  spectral  notching  is  added  during  generation  of  the  communication  sym¬ 
bols,  there  is  a  dramatic  improvement  in  Pe  and  Pb  performance.  Figure  4.11 
and  Fig.  4.12  show  Pe  and  Pb  versus  normalized  signal-to-noise  ratio,  respectively, 
for  this  scenario.  System  performance  when  spectral  notching  is  employed  nearly 
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Figure  4.10  Jamming  Present  with  No  Spectral  Notching :  Pb  vs  E^/Nq  for  TD- 
MPSK  Signaling  with  J/S  =  3.14  dB  and  P  =  31  Sinusoids 

matches  system  performance  when  no  jamming  is  present  (compare  these  two  fig¬ 
ures  with  Fig.  4.5  and  Fig.  4.6).  However,  because  the  jamming  power  is  not  entirely 
removed,  there  is  a  slight  increase  in  error  by  comparison  with  scenarios  containing 
no  jamming  (see  TD-BPSK  and  TD-QPSK  performance  at  higher  normalized  signal- 
to-noise  ratios).  At  higher  normalized  signal-to-noise  ratios  there  is  a  separation  in 
the  error  between  the  two  modulation  types.  If  no  jamming  had  been  present,  Pb  for 
the  two  modulation  types  would  have  been  identical.  However,  since  the  signal-to- 
interference  ratio  is  held  constant,  the  interference  affects  higher  modulation  levels, 
such  as  TD-QPSK,  more  than  lower  levels  of  modulation,  such  as  TD-BPSK. 


4.7  Multiple  Access  Interference:  Njj  =  3  Users  With  No  Spectral 
Notching 

The  next  scenario  involved  multiple  access  interference  given  only  one  addi¬ 
tional  user  (transmitter)  is  present.  This  additional  transmitter  uses  the  same  mag¬ 
nitude  frequency  spectrum  as  the  transmitter /receiver  pair  (link)  of  interest.  The 
multiple  access  phase  coding  of  the  additional  transmitter  is  random  relative  to  the 
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E/ryffi) 

Figure  4.11  Jamming  Present  With  Spectral  Notching:  Pe  vs  E^/Nq  for  TD-MPSK 
Signaling  with  J/S  =  3.14  dB  and  P  =  31  Sinusoids 


WdB) 


Figure  4.12  Jamming  Present  With  Spectral  Notching:  Pb  vs  E^/Nq  for  TD-MPSK 
Signaling  with  J/S  =  3.14  dB  and  P  =  31  Sinusoids 

multiple  access  phase  coding  used  for  the  link  of  interest.  Figure  4.13  and  Fig.  4.14 
show  Pe  and  Pb  versus  normalized  signal-to-noise  ratio,  respectively,  for  this  sce¬ 
nario.  In  this  case,  interference  from  the  additional  transmitter  can  be  modeled 
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as  being  broad  band  and  has  the  same  power  as  the  link  of  interest.  The  signal- 
to-interference  ratio  remains  constant  for  different  levels  of  modulation.  As  in  the 
the  narrow  band  jamming  scenario,  the  interference  affects  the  higher  modulation 
levels  more  than  the  lower  modulation  levels.  Once  again  this  is  most  obvious  by 
comparing  Pb  curves  for  the  TD-QPSK  and  TD-BPSK  data  modulations  where  the 
multiple  access  interference  affects  TD-QPSK  more  than  TD-BPSK. 


Figure  4.13  Multiple  Access  Interference  for  Nu  =  3  Users  With  No  Spectral 
Notching:  Pe  vs  Eb/N0  for  TD-MPSK  Signaling  with  P  =  31  Si¬ 
nusoids 


4.8  Multiple  Access  Interference:  N u  =  2  to  32  Users  With  No  Spectral 
Notching 

As  the  number  of  users  (Ny) increases  the  probability  of  error  increases.  Inter¬ 
ference  by  Nu  users  at  a  particular  Eb/Ni  is  calculated  using  (4.3),  where  {Nv  —  2) 
is  the  number  of  additional  users,  and  it  is  assumed  that  the  spectrum  magnitudes 
of  all  users  are  identical.  Equation  (4.3)  is  a  modified  version  of  (3.27  which  only 
took  into  consideration  the  effect  of  one  additional  user. 
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Figure  4.14  Multiple  Access  Interference  for  Njj  =  3  Users  With  No  Spectral 
Notching :  Pb  vs  Eb/N0  for  TD-MPSK  Signaling  with  P  =  31  Si¬ 
nusoids 


Ni  2R(Nu-  2)^|(a(„)s„)2| 

The  result  of  (4.3)  is  substituted  into  the  equation  for  calculating  Eb/(N0  +  Nj) 
which  in  turn  is  used  in  the  analytic  Pb  and  Pe  equations.  As  Njj  increases  the 
probability  of  error  increases.  Figure  4.15  and  Fig.  4.16  show  Pe  and  Pb  versus 
Njj,  respectively,  for  Ay-  =  2  to  32  users  and  TD-QPSK  data  modulation.  For 
these  simulations,  the  normalized  signal-to-noise  ratio  was  maintained  constant  at 
Eb/N0  =  6.0  dB.  The  interference  caused  by  the  additional  users  is  proportional  to 
the  number  of  additional  users  present  in  the  network. 

4.9  Asynchronous  vs.  Synchronous  Multiple  Access  Performance 

For  equal  amplitude  sinusoids,  asynchronous  cross-correlation  is  approximately 
the  same  as  synchronous  cross-correlation.  Figure  4.17  and  Fig.  4.18  show  Pe  and 
Pb  versus  Njj  respectively  for  asynchronous  networks.  A  comparison  with  Fig.  4.15 
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Figure  4.15  Multiple  Access  Interference  for  NE  =  2  to  32  Users  With  No  Spectral 
Notching :  PE  vs  Nrj  for  TD-QPSK  Signaling  with  P  =  31  Sinusoids 


Figure  4.16  Multiple  Access  Interference  for  NE  =  2  to  32  Users  With  No  Spectral 
Notching:  PE  vs  NE  for  TD-QPSK  Signaling  with  P  =  31  Sinusoids 

and  Fig.  4.16  shows  bit  error  performance  is  nearly  identical  under  multiple  access 
conditions  in  synchronous  and  asynchronous  networks. 
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Figure  4.17  Multiple  Access  Interference  with  Njj  =  2  to  32  Asynchronous  Users: 
PE  vs  Njj  for  TD-QPSK  Signaling  with  P  =  31  Sinusoids 


Figure  4.18  Multiple  Access  Interference  with  Njj  =  2  to  32  Asynchronous  Users: 
PB  vs  Njj  for  TD-QPSK  Signaling  with  P  =  31  Sinusoids 

4.10  Multiple  Access  Interference:  Njj  =  3  Users  With  Spectral  Notch¬ 
ing 

When  an  additional  multiple  access  user  is  added  to  a  network  employing 
spectral  notching,  performance  is  degraded.  As  simulated  here,  additional  multiple 
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access  users  are  assumed  to  be  experiencing  the  same  electromagnetic  interference 
(spectral  location  and  power)  as  the  link  of  interest.  Therefore,  all  network  users 
generate  identical  spectral  notches  as  the  link  of  interest.  Figure  4.19  and  Fig.  4.20 
show  Pb  and  Pe  versus  normalized  signal-to-noise  ratio,  respectively,  for  the  Njj  =  3 
user  case.  As  indicated,  the  performance  of  a  network  using  notching  is  degraded 
when  an  additional  user  is  added.  Furthermore,  the  degradation  indicated  is  greater 
than  observed  in  a  system  without  notching  (see  Fig.  4.13  and  Fig.  4.14  for  compar¬ 
ison).  This  occurs  because  the  notched  signal  power  is  spread  over  fewer  frequency 
components  than  the  signal  without  notching;  fewer  frequency  components  yields 
higher  cross-correlation  with  additional  users  and  produces  higher  bit  error  rates. 


Figure  4.19  Multiple  Access  Interference  for  Nu  =  3  Users  With  Spectral  Notching-. 
PE  vs  Eb/N0  for  TD-MPSK  Signaling  with  P  =  31  Sinusoids 


4.11  Multiple  Access  Interference:  N u  =  2  to  32  Users  With  Spectral 
Notching 

As  in  the  multiple  access  case  with  no  spectral  notching,  the  amount  of  multiple 
access  interference  is  proportional  to  the  number  of  additional  users.  However,  the 
amount  of  interference  caused  by  each  user  is  greater  than  experienced  in  the  case 
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Figure  4.20  Multiple  Access  Interference  for  Nu  =  3  Users  With  Spectral  Notching-. 
PB  vs  Eb/N0  for  TD-MPSK  Signaling  with  P  =  31  Sinusoids 


with  no  spectral  notching.  Figure  4.21  and  Fig.  4.22  show  PB  and  PB  versus  Njj, 
respectively,  for  NB  =  2  to  32  users  with  spectral  notching.  A  comparison  with 
results  presented  in  Fig.  4.15  and  Fig.  4.16  show  how  probability  of  error  increases 
at  a  greater  rate  for  increasing  Nu  when  notching  is  present. 

As  the  spectral  notch  width  increases,  multiple  access  interference  increases. 
Figure  4.23  and  Fig.  4.24  show  PB  and  PB  versus  notch  width  Ns  (the  number  of 
sinusoidal  components  zeroed),  respectively,  for  Nu  =  6  users  and  TD-QPSK  data 
modulation.  The  multiple  access  interference  for  these  simulations  is  from  4  addi¬ 
tional  users  having  identical  magnitude  spectra  at  the  transmitters.  The  simulation 
process  started  by  modeling  all  transmitters  as  using  P  =  31  sinusoids;  PB  and  PB 
performance  were  then  established  with  the  multiple  access  interference  present.  The 
highest  frequency  sinusoid  was  then  removed  and  PB  and  PB  again  determined  with 
the  multiple  access  interference  present.  The  process  was  repeated  until  only  the 
lowest  frequency  sinusoid  remained.  Clearly,  as  the  notch  width  increases  both  PB 
and  PB  increase. 
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Figure  4.21  Multiple  Access  Interference  for  Njj  =  2  to  32  Users  With  Spectral 
Notching :  Synchronous  Users:  Pe  vs  Njj  for  TD-QPSK  Signaling 
with  P  =  31  Sinusoids 


Figure  4.22  Multiple  Access  Interference  for  Njj  =  2  to  32  Users  With  Spectral 
Notching:  Synchronous  Users:  Pb  vs  Njj  for  TD-QPSK  Signaling 
with  P  =  31  Sinusoids 

4.12  Jamming  and  Nu  —  3  Users  Present  With  No  Spectral  Notching 

When  jamming  and  Njj  =  3  Users  are  present,  the  interference  introduced 
into  the  system  effectively  compounds.  Figure  4.25  and  Fig.  4.26  show  P&  and 
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Figure  4.23  Multiple  Access  Interference  for  NB  =  6  Users  With  Spectral  Notching : 
PE  vs  Notch  Width  Ns  for  TD-QPSK  Signaling 


Figure  4.24  Multiple  Access  Interference  for  Njj  =  6  Users  With  Spectral  Notching : 
PB  vs  Notch  Width  Ns  for  TD-QPSK  Signaling 

Pb  versus  the  normalized  signal-to-noise  ratio,  respectively,  for  the  case  when  both 
jamming  and  Nv  =  3  users  are  present.  As  expected,  the  degradation  induced  in 


4-18 


system  performance  is  greater  than  what  occurs  when  either  the  additional  user  or 
the  jammer  are  independently  introduced. 


E/ryffi) 

Figure  4.25  Jamming  and  TVjy  =  3  Users  Present  With  No  Spectral  Notching :  Pe 
vs  Eb/N0  for  TD-MPSK  Signaling  with  .J/S  =  3.14  dB  and  P  =  31 
Sinusoids 


WdB) 

Figure  4.26  Jamming  and  Nv  =  3  Users  Present  With  No  Spectral  Notching :  Pb 
vs  Eb/N0  for  TD-MPSK  Signaling  with  ,J/S  =  3.14  dB  and  P  =  31 
Sinusoids 
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4.13  Jamming  and  N jj  2  to  32  Users  Present  With  No  Spectral  Notching 

When  jamming  is  present  it  has  the  affect  of  raising  the  noise  floor  in  a  mul¬ 
tiple  access  environment.  Figure  4.27  and  Fig.  4.28  show  Pe  and  Pb  versus  Njj, 
respectively,  for  J/S  =  3.14  dB  and  TD-QPSK  data  modulation.  By  comparison 
with  Fig.  4.15  and  Fig.  4.16,  Pe  and  Pb  are  greater  than  the  Njj  =  3  multiple  access 
case  with  no  jamming  present. 


Figure  4.27  Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  32  Users 
No  Spectral  Notching :  Synchronous  Users:  Pe  vs  Njj  for  TD-QPSK 
Signaling  with  P  =  31  Sinusoids 

4.14  Jamming  and  Njj  —  3  Users  Present  With  Spectral  Notching 

When  jamming,  an  additional  user,  and  a  notch  are  all  present,  the  probability 
of  symbol  error  and  probability  of  bit  error  performance  are  almost  as  good  as  when 
an  additional  user  is  present  and  a  notch  is  present.  Figure  4.29  and  Fig.  4.30 
show  Pe  and  Pb  versus  normalized  signal  to  noise  ratio  (Eb/N0),  respectively.  As 
indicated,  spectral  notching  is  able  to  remove  most  of  the  jamming,  yet  there  remains 
some  degradation  caused  by  the  notch.  When  adding  a  notch  to  remove  jamming 
the  multiple  access  interference  should  be  considered.  It  is  possible  to  remove  most 
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Figure  4.28  Jamming  and  Multiple  Access  Interference  for  Nu  =  2  to  32  Users 
No  Spectral  Notching-.  Synchronous  Users:  Pg  vs  Njj  for  TD-QPSK 
Signaling  with  P  =  31  Sinusoids 

of  the  effects  of  jamming  but  increase  the  effect  of  multiple  access  interference  to 
such  a  degree  that  that  the  system  is  worse  off.  A  notch  should  be  large  enough  to 
remove  most  of  the  jamming.  However,  it  should  not  be  so  large  that  it  increases 
the  multiple  access  interference  to  such  a  degree  that  the  system  is  worse  off. 

4.15  Jamming  and  Nv  —  2  to  32  Users  Present  With  With  Spectral 
Notching 

Jamming  has  the  affect  of  raising  the  noise  floor  in  a  multiple  access  environ¬ 
ment.  Figure  4.31  and  Fig.  4.32  show  Pe  and  Pb  versus  Nu,  respectively,  for  the 
number  of  users.  The  performance  of  this  scenario  is  degraded  slightly  versus  the 
scenario  with  a  notch  but  no  jamming.  This  is  because  there  is  still  some  residual 
jamming  that  wasn’t  notched  away.  It  is  interesting  to  compare  this  scenario  versus 
the  scenario  with  jamming  but  no  notch.  Figure  4.33  shows  this  comparison.  After 
five  additional  users,  probability  of  bit  error  is  better  for  the  case  without  the  notch. 
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Figure  4.29  Jamming  and  Nu  =  3  Users  Present  With  With  Spectral  Notching :  Pb 
vs  Eb/N0  for  TD-MPSK  Signaling  with  J/S  =  3.14  dB  and  P  =  31 
Sinusoids 


WdB) 


Figure  4.30  Jamming  and  Nu  =  3  Users  Present  With  With  Spectral  Notching-.  Pb 
vs  Eb/N0  for  TD-MPSK  Signaling  with  ,J/S  =  3.14  dB  and  P  =  31 
Sinusoids 


4.16  Summary 

Simulated  results  for  probability  of  symbol  error  (Pe)  and  probability  of  bit 
error  (Pb)  performance  versus  normalized  signal-to-noise  ratio  ( Eb/N0 )  are  shown 
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Figure  4.31  Jamming  and  Multiple  Access  Interference  for  Nu  =  2  to  32  Users 
With  Spectral  Notching :  Synchronous  Users:  Pe  vs  Nv  for  TD-QPSK 
Signaling  with  P  =  31  Sinusoids 


Figure  4.32  Jamming  and  Multiple  Access  Interference  for  Nu  =  2  to  32  Users 
With  Spectral  Notching-.  Synchronous  Users:  Pb  vs  Nu  for  TD-QPSK 
Signaling  with  P  =  31  Sinusoids 

consistent  (nearly  identical)  to  analytic  results  obtained  from  equations  provided  in 
Section  2.2.1  and  Section  2.2.2.  This  includes  the  cases  when  only  spectral  notching 
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X  Analytic  QPSK  (Jamming) 

•  Analytic  QPSK  (Jamming  and  Notching) 

10  4I -  1 

2  7  12  17  22  27  32 

Number  of  Users  (N  ) 

Figure  4.33  Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  32  Users 
With  and  without  Spectral  Notching :  Pb  vs  Nu  for  TD-MPSK  Signal¬ 
ing  with  P  =  31  Sinusoids 

or  only  jamming  are  considered,  as  well  as,  the  case  when  jamming  with  spectral 
notching  is  considered.  Results  for  scenarios  containing  multiple  access  interference 
are  consistent  with  analytic  expressions  as  well.  Again,  this  includes  multiple  access 
interference  scenarios  where  only  spectral  notching  or  only  jamming  are  present,  as 
well  as,  the  multiple  access  case  when  jamming  with  spectral  notching  is  considered. 
Figure  4.34  shows  for  an  Eb/N0  of  6  dB,  an  equal  power  additional  user  degrades  Pb 
from  7.9  x  10”2  to  8.5  x  1CU2.  If  an  equal  power  additional  user  and  jamming  are 
added  to  the  system,  Pb  further  degrades  to  9.8  x  1CT2.  If  a  spectral  notch  is  added  to 
the  system  when  jamming  is  not  present,  Pb  degrades  to  8.8  x  10~2.  This  is  because 
the  communication  symbols  are  not  spread  over  as  many  spectral  components,  which 
increases  the  correlation  between  different  users,  which  then  increases  the  multiple 
access  interference.  If  jamming,  an  additional  equal  power,  user  and  a  spectral  notch 
are  present,  Pb  is  8.9  x  1CU2.  Pb  in  the  case  with  an  additional  user,  jamming,  and 
spectral  notching  is  slightly  larger  than  the  case  with  an  additional  user  and  notching 
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because  the  small  amount  of  jamming  (6  %)  that  was  not  notched  away  causes  some 
degradation. 


Figure  4.34 


Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  3  Users  With 
and  without  Spectral  Notching-.  Pb  vs  E^/Nq  for  TD-QPSK  Signaling 
with  P  =  31  Sinusoids 
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5.  Modeling  and  Simulation:  Orthogonal  Code  Selection 

5.1  Introduction 

This  chapter  provides  an  overview  of  modeling  parameters  and  scenarios  used 
for  testing  performance  of  orthogonally  coded  users  in  the  proposed  TD-MPSK  sys¬ 
tem.  The  parameters  described  include  the  communication  symbol  spectral  response, 
the  jammer  spectral  response  and  the  notched  symbol  spectral  response.  The  mod¬ 
eled  scenarios  include:  1)  multiple  access  interference,  2)  multiple  access  interference 
with  spectral  notching,  3)  multiple  access  and  jamming  interference  without  spec¬ 
tral  notching  and  4)  multiple  access  and  jamming  interference  with  spectral  notching. 
The  scenarios  considered  help  characterize  the  TD-MPSK  system’s  ability  to  commu¬ 
nicate  in  the  presence  of  noise  and  interference  when  the  system  is  using  orthogonally 
coded  users. 

5.2  Modeling  Parameters 

5.2.1  Spectral  Response  of  Communication  Symbol.  For  the  or¬ 
thogonal  code  selection  case,  the  spectral  response  is  almost  the  same  as  presented 
in  the  previous  chapter  in  Fig.  4.1.  The  only  difference  is  the  multiple  access  phase 
codes  of  the  orthogonal  users  are  interrelated.  Equation  (5.1)  is  a  discrete  form 
of  the  analytic  representation  of  communication  symbols  in  an  orthogonal  network. 
Equation  (5.2)  describes  the  phase  codes  for  the  different  users  in  the  network.  Equa¬ 
tions  (5.3)  and  (5.4)  describe  the  variable  G  and  the  number  of  possible  users  in  the 
network.  Each  member  in  the  set  of  orthogonal  users  uses  the  same  <fip  to  adjust 
the  phase  of  each  of  the  sinusoids.  The  phase  of  each  sinusoid  for  each  symbol  is 
further  adjusted  based  on  user  number  (va)  where  va  =  0  identifies  the  primary  user 
being  modeled.  The  remaining  orthogonal  users  are  identified  as  additional  users. 
The  first  additional  orthogonal  user  is  based  on  using  v0  =  1,  the  second  is  based 
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on  using  v0  =  2,  and  so  on,  up  to  the  thirtieth  additional  user  which  is  based  on 
v0  =  30. 

4^  ^  =  h  Ap  cos  +  ^  +  0kV))  v  =  0,...,G-l  (5.1) 

p=  i 

0p]  =  4>p  +  s§n(A)  ^4  (5-2) 

9=1 

31 

G  =  y^sgn(Ar)  (5.3) 

r= 1 

l  t  >  o 

sgn(t)  =  <  0  t  =  0  (5.4) 

-1  t  <  0 

5.2.2  Spectral  Response  of  Jamming.  The  spectrum  of  the  jamming 

is  the  same  as  that  of  the  previous  chapter. 

5.2.3  Spectrum  of  Notched  User.  The  magnitude  of  the  spectrum  of 
the  notched  users  is  the  same  as  that  of  the  previous  chapter. 

5.3  Synchronous  Multiple  Access  Interference:  Orthogonal  Users  With 
No  Spectral  Notching 

Probability  of  error  (Pe  and  Pb)  generally  increases  as  the  number  of  multiple 
access  users  (Njj)  increases.  However,  for  the  case  of  a  synchronous  network  con¬ 
taining  orthogonal  users,  Pe  and  Pb  are  unaffected  as  Nv  increases  provided  each 
additional  user  is  mutually  orthogonal  to  all  previous  users.  Figure  5.1  and  Fig.  5.2 
show  Pe  and  Pb  versus  Nu,  respectively,  for  Ay/  =  2  to  32  users  and  TD-QPSK 
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data  modulation.  For  these  simulations,  the  normalized  signal-to-noise  ratio  was 
maintained  constant  at  E^/Nq  =  6.0  dB.  As  indicated  by  the  data  in  these  figures, 
for  P  =  31  sinusoids  as  given  in  (5.1),  it  is  possible  to  assign  multiple  access  phase 
coding  in  accordance  with  5.2  through  5.4  such  that  G  users  in  the  network  are 
mutually  orthogonal. 


2  7  12  17  22  27  32 


Number  of  Users  (N  ) 

Figure  5.1  Multiple  Access  Interference  for  Nv  =  2  to  32  Orthogonal  Users  With 
No  Spectral  Notching :  Pe  vs  Nv  for  TD-QPSK  Signaling  with  P  =  31 
Sinusoids 

5.4  Asynchronous  Multiple  Access  Interference:  Orthogonal  Users  With 
No  Spectral  Notching 

The  cross-correlation  of  orthogonal  users  is  approximately  the  same  as  asyn¬ 
chronous  cross-correlation  with  randomly  coded  users  with  the  same  magnitude  spec¬ 
trum  as  the  orthogonal  users.  Figure  5.3  and  Fig.  5.4  show  Pe  and  Pb  versus  Njj 
respectively  for  asynchronous  networks. 
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Figure  5.2  Multiple  Access  Interference  for  Njj  =  2  to  32  Orthogonal  Users  With 
No  Spectral  Notching:  PB  vs  Nrj  for  TD-QPSK  Signaling  with  P  =  31 
Sinusoids 
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Figure  5.3  Multiple  Access  Interference  with  Njj  =  2  to  32  Asynchronous  Users: 
Pe  vs  Njj  for  TD-QPSK  Signaling  with  P  =  31  Sinusoids 

5.5  Multiple  Access  Interference:  NB  =  2  to  32  Users  With  Spectral 
Notching 

If  spectral  notching  is  introduced  into  a  group  of  orthogonal  users,  orthogo¬ 
nality  is  lost.  Figure  5.5  and  Fig.  5.6  show  Pe  and  PB  versus  the  number  of  users, 
respectively,  when  spectral  notching  is  efhjiloyed. 


Figure  5.4  Multiple  Access  Interference  with  NB  =  2  to  32  Asynchronous  Users: 
PB  vs  Njj  for  TD-QPSK  Signaling  with  P  =  31  Sinusoids 


Figure  5.5  Multiple  Access  Interference  for  NB  =  2  to  32  Orthogonal  Users  With 
Spectral  Notching :  PB  vs  NB  for  TD-QPSK  Signaling  with  P  =  31 
Sinusoids 

As  the  spectral  notch  width  increases,  multiple  access  interference  increases. 
Figure  5.7  and  Fig.  5.8  show  PB  and  PB  the  versus  spectral  notch  width,  respectively. 
The  multiple  access  interference  for  the  simulations  comes  from  four  additional  or- 
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-  Analytic  QPSK 

•  TD-QPSK 
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Figure  5.6  Multiple  Access  Interference  for  Njj  =  2  to  32  Orthogonal  Users  With 
Spectral  Notching :  PB  vs  Nv  for  TD-QPSK  Signaling  with  P  =  31 
Sinusoids 


thogonal  users.  Each  of  these  orthogonal  users  has  the  same  magnitude  spectrum 
as  the  transmitter /receiver  pair.  The  simulations  start  by  modeling  performance 
when  the  transmitter  uses  all  31  sinusoids.  At  this  point,  users  maintain  mutual 
orthogonality  and  do  not  induce  any  multiple  access  interference  during  the  correla¬ 
tion  process;  any  resultant  symbol  and  bit  errors  are  due  solely  to  the  environmental 
noise.  At  the  next  step,  only  the  highest  frequency  sinusoid  is  removed  and  PE  and 
Pb  determined  again.  At  this  point  the  orthogonal  users  lost  their  orthogonality. 
The  process  is  repeated  by  removing  one  sinusoid  at  a  time  until  there  is  only  one  si¬ 
nusoid  remaining.  As  the  spectral  notch  width  increases  (by  progressively  removing 
one  sinusoid  at  a  time)  PE  and  PB  increase. 


5.6  Jamming  and  Orthogonal  Users  Present  With  No  Spectral  Notching 

As  in  the  randomly  coded  user  case,  jamming  has  a  tendency  to  raise  the 
noise  floor.  Figure  5.9  and  Fig.  5.10  show  PE  and  PB  versus  NE,  respectively,  when 
jamming  and  orthogonal  users  are  present  and  no  spectral  notching  is  applied.  For 
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Figure  5.7  Multiple  Access  Interference  for  Njj  =  6  Orthogonal  Users  With  Spectral 
Notching :  PE  vs  Spectral  Notch  Width  Ns  for  TD-QPSK  Signaling 


Figure  5.8  Multiple  Access  Interference  for  Njj  =  6  Orthogonal  Users  With  Spectral 
Notching :  PE  vs  Spectral  Notch  Width  Ns  for  TD-QPSK  Signaling 

the  Nu  =  2  to  32  case,  PE  and  PE  remain  constant.  As  expected,  the  probability  of 
error  is  greater  here  than  experienced  in  the  case  when  no  jamming  present. 
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Figure  5.9  Jamming  and  Multiple  Access  Interference  for  Nu  =  2  to  32  Orthogonal 
Users  No  Spectral  Notching :  PE  vs  Nv  for  TD-QPSK  Signaling  with 
P  =  31  Sinusoids 


Figure  5.10  Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  32  Orthog¬ 
onal  Users  No  Spectral  Notching-.  PE  vs  NE  for  TD-QPSK  Signaling 
with  P  =  31  Sinusoids 


5.7  Jamming  and  Orthogonal  Users  Present  With  Spectral  Notching 


Although  adding  a  spectral  notch  to  the  communication  system  removes  most 
of  the  jamming,  it  degrades  and  possibly  destroys  the  desired  mutual  orthogonality 
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between  users.  Figure  5.11  and  Fig.  5.12  show  Pe  and  Pb  versus  the  number  of 
users,  respectively,  when  jamming  and  orthogonal  users  are  present  and  spectral 
notching  is  applied.  Results  in  these  figures  show  slight  degradation  in  performance 
by  comparison  with  the  case  where  a  spectral  notch  was  used  without  jamming 
present.  This  is  because  there  is  some  residual  jamming. 

There  is  a  trade-off  that  needs  to  be  made  between  spectral  notching  and 
maintaining  orthogonality.  If  there  are  only  few  additional  orthogonal  users,  spectral 
notching  might  be  of  some  benefit.  Although  orthogonality  is  lost  in  this  case, 
the  jamming  power  that  is  removed  as  a  result  of  the  spectral  notch  may  make 
this  trade-off  a  viable  option.  In  cases  where  there  are  a  relatively  high  number  of 
orthogonal  users,  multiple  access  interference  that  remains  after  spectral  notching 
becomes  greater  than  the  jamming  interference  that  is  actually  removed  by  the 
spectral  notch.  For  the  specific  scenarios  modeled  in  this  chapter,  once  the  number  of 
orthogonal  users  exceeds  two,  notching  causes  more  degradation  than  jamming.  The 
loss  of  orthogonality  caused  by  the  notching  causes  the  multiple  access  interference 
to  exceed  the  jamming  interference. 

5.8  Special  Case:  Reassigning  Phase  Codes  to  Maintain  Orthogonality 

Spectral  notching  removes  some  of  the  sinusoids  making  up  the  symbol.  How¬ 
ever,  if  the  number  of  users  is  reduced  to  the  number  of  available  sinusoids,  the 
multiple  access  code  phases  can  be  reassigned  to  restore  orthogonality.  The  new 
multiple  access  code  phases  are  calculated  using  (5.1)  through  (5.4).  For  the  case 
of  spectral  notching  discussed  in  Section  4.2.3,  the  value  of  G  drops  from  31  to  23 
because  of  the  number  of  frequency  components  that  were  notched  away. 

Figure  5.13  and  Fig.  5.14  show  Pe  and  Pb  versus  the  number  of  users,  respec¬ 
tively,  when  coding  phases  have  been  reassigned  and  spectral  notching  is  applied. 
The  number  of  remaining  sinusoids  that  make  up  a  symbol  is  23  which  allows  the 
creation  of  23  orthogonal  users. 
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Figure  5.11  Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  32  Orthogo¬ 
nal  Users  With  Spectral  Notching:  Pe  vs  Njj  for  TD-QPSK  Signaling 
with  P  =  31  Sinusoids 


Figure  5.12  Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  32  Orthogo¬ 
nal  Users  With  Spectral  Notching:  Pb  vs  Njj  for  TD-QPSK  Signaling 
with  P  =  31  Sinusoids 

Figure  5.15  and  Fig.  5.16  show  PE  and  PE  versus  the  number  of  users,  respec¬ 
tively,  when  coding  phases  have  been  reassigned,  jamming  is  present,  and  spectral 
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Figure  5.13  Multiple  Access  Interference  for  Ny  =  2  to  24  Orthogonal  Users  With 
Spectral  Notching,  Orthogonality  is  Restored :  Pe  vs  Nv  for  TD-QPSK 
Signaling  with  P  =  31  Sinusoids 
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Figure  5.14  Multiple  Access  Interference  for  Njj  =  2  to  24  Orthogonal  Users  With 
Spectral  Notching,  Orthogonality  is  Restored:  Pb  vs  Njj  for  TD-QPSK 
Signaling  with  P  =  31  Sinusoids 

notching  is  applied.  As  in  the  previous  case  there  are  23  orthogonal  users.  Perfor¬ 
mance  versus  the  previous  case  is  degraded  since  some  of  the  jamming  remains  in 
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the  frequency  components  that  were  not  notched.  However,  this  case  is  only  slightly 
degraded  versus  the  previous  case  because  the  majority  of  the  jamming  was  notched 
away. 
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Figure  5.15  Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  24  Orthogo¬ 
nal  Users  With  Spectral  Notching,  Orthogonality  Restored :  Pe  vs  Njj 
for  TD-QPSK  Signaling  with  P  =  31  Sinusoids 


5.9  Summary 

The  cross-correlation  between  communication  symbols  of  different  synchronous 
users  can  be  made  identically  zero  through  proper  selection  of  multiple  access  phase 
codes  (orthogonal  signaling).  For  a  synchronous  network  containing  orthogonal 
users,  symbol  and  bit  error  probabilities  are  unaffected  as  the  number  of  orthog¬ 
onal  network  users  increases.  The  addition  of  jamming  does  not  change  the  amount 
of  multiple  access  interference  although  it  does  increase  the  probability  of  error.  The 
introduction  of  spectral  notching  increases  the  multiple  access  interference  since  net¬ 
work  users  in  the  network  are  no  longer  orthogonal.  If  the  multiple  access  phase 
codes  are  properly  reassigned  after  spectral  notching,  a  reduced  number  of  orthog¬ 
onal  users  can  be  created.  Figure  5.17  summarizes  the  affect  of  jamming,  spectral 
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Figure  5.16  Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  24  Orthogo¬ 
nal  Users  With  Spectral  Notching,  Orthogonality  Restored :  Pb  vs  Njj 
for  TD-QPSK  Signaling  with  P  =  31  Sinusoids 

notching,  and  reassigning  phase  codes  to  restore  orthogonality.  When  no  jamming  is 
present,  Pb  is  2.6  x  10-4.  When  jamming  is  added  to  the  system,  Pb  error  increases 
sharply  to  5.1  x  10~3.  When  a  spectral  notch  is  added  to  the  system  to  remove 
the  jamming,  Pb  increases  to  approximately  2.5  x  10“2  for  Njj  =  4  to  30.  The 
notch  removed  most  of  the  jamming  interference  but  increased  the  multiple  access 
interference.  When  the  spectral  codes  were  reassigned,  Pb  decreased  dramatically 
to  4.0  x  10”4.  The  slight  increase  in  error  versus  the  case  with  no  jamming  and  no 
spectral  notching  is  due  to  the  small  amount  of  jamming  (6%)  that  was  not  notched 
away. 
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Figure  5.17 


Jamming  and  Multiple  Access  Interference  for  Njj  =  2  to  32  Orthog¬ 
onal  Users  With  and  without  Spectral  Notching :  Pb  vs  Njj  for  TD- 
QPSK  Signaling  with  P  =  31  Sinusoids 
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6.  Conclusions  and  Recommendations 


6.1  Summary 

First  and  foremost,  this  thesis  provides  the  introduction,  development  and 
characterization  of  a  spectrally  encoded  transform  domain,  M-Ary  phase  shift  keying 
(TD-MPSK)  technique  that  provides  both  multiple  access  capability  and  interference 
suppression  (avoidance).  As  presented,  the  proposed  TD-MPSK  system  uses  a  form 
of  spectral  encoding,  i.e.,  the  application  of  independent  data  and  multiple  access 
phase  modulations  to  sinusoidal  spectral  components,  to  generate  multiple  access 
communication  symbols  which  are  subsequently  demodulated  in  the  time  domain 
using  a  conventional  correlation  receiver.  The  previous  transform  domain  communi¬ 
cation  systems  considered  relied  primarily  on  binary  phase  shift  keying  (BPSK)  and 
cyclic  shift  keying  (CSK)  to  provide  data  modulation.  This  thesis  provides  a  review 
of  conventional  MPSK  signaling,  transform  domain  communications,  and  code  divi¬ 
sion  multiple  access.  This  background  is  followed  by  a  description  of  the  TD-MPSK 
symbol  definition  and  phase  coding  used  for  data  modulation  and  multiple  access  im¬ 
plementation.  Performance  is  characterized  in  terms  of  multiple  access  interference 
and  jamming  interference.  The  probability  of  error  (symbol  and  bit)  for  a  spectrally 
encoded  TD-MPSK  system  is  discussed  in  terms  of  different  types  of  interference 
present.  Monte  Carlo  simulation  techniques  are  used  to  verify  the  analytical  expres¬ 
sions  derived  explicitly  for  the  proposed  TD-MPSK  system.  Various  modeling  and 
simulation  scenarios  are  considered,  including,  those  containing  jamming  interfer¬ 
ence,  multiple  access  interference,  and  combinations  thereof,  both  with  and  without 
TD-MPSK  spectral  notching  employed.  This  thesis  concludes  with  a  description  of 
Monte  Carlo  simulations  conducted  to  verify  TD-MPSK  multiple  access  performance 
with  orthogonal  users. 
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6.2  Conclusions 


6.2.1  Performance  in  the  Presence  of  Interference.  This  thesis 
shows  that  probability  of  symbol  error  (Pe)  and  probability  of  bit  error  (Pb)  versus 
normalized  signal-to-noise  ratio  for  the  proposed  TD-MPSK  system  can  be  reliably 
estimated  using  conventional  MPSK  error  equations  for  an  orthonormal  signal  space, 
i.e.,  a  pair  of  unit  energy  orthogonal  functions  that  mathematically  span  the  two- 
dimensional  signal  space.  Both  Pe  and  Pb  can  be  reliably  predicted  in  the  presence 
of  additional  network  users  (multiple  access  interference),  intentional  jamming  and 
spectral  notching.  The  error  estimates  derived  for  different  scenarios  were  validated 
using  Monte  Carlo  simulation  and  analysis. 

6.2.2  Multiple  Access  Performance  (Randomly  Coded  Users). 

The  amount  of  multiple  access  interference  induced  by  randomly  coded  network 
users  can  be  estimated  using  equations  derived  as  part  of  this  research.  Both  Pe 
and  Pb  versus  the  number  of  users  can  be  predicted  using  the  derived  multiple  access 
interference  expressions  and  previously  derived  conventional  MPSK  error  equations. 
Interference  effects  on  Pe  and  Pb  can  also  be  predicted  in  the  presence  of  jamming 
and  spectral  notching.  Error  estimates  were  validated  for  different  scenarios  using 
Monte  Carlo  simulation  and  analysis. 

6.2.3  Multiple  Access  Performance  (Orthogonally  Coded  Users). 

Multiple  access  interference  in  a  synchronous  network  can  be  eliminated  if  users  are 
made  mutually  orthogonal.  For  a  limited  number  of  users,  the  TD-MPSK  multiple 
access  phase  codes  and  sinusoidal  component  amplitudes  can  be  selected  to  provide 
this  orthogonality.  If  the  orthogonally  coded  users  are  used  in  an  asynchronous 
network,  error  performance  degrades  and  approaches  that  of  randomly  coded  users 
operating  in  an  asynchronous  network. 


6-2 


6.3  Recommendations  for  Future  Research 


6.3.1  Quadrature  Amplitude  Modulation.  The  TD-MPSK  symbols 
developed  as  part  of  this  research  are  fundamentally  based  on  two  orthogonal  signals 
defining  (mathematically  spanning)  the  two-dimensional  signal  space.  By  combining 
amplitude  modulation  of  the  orthogonal  symbols  (?/q  and  1^2)  with  phase  shift  keying 
(PSK),  a  quadrature  amplitude  modulation  (QAM)  system  can  be  developed. 

6.3.2  Synchronization.  Research  has  previously  been  conducted  on 
synchronization  of  a  TDCS  system  using  BPSK  and  CSK  modulations  [2],  This 
research  could  be  reviewed  to  see  how  it  applies  to  a  TD-MPSK  system. 

6.3.3  Orthogonal  Jammers.  One  possible  use  of  this  research  is  the 
creation  of  a  combination  jammer/communication  system.  This  combination  system 
could  generate  a  jammer  waveform  which  is  orthogonal  (non-interfering)  to  friendly 
communication  systems  but  disruptive  to  enemy  communication  and  radar  systems. 

6.3.4  Propagation  Effects.  Further  research  can  be  done  on  propa¬ 
gation  effects  on  a  TDCS  system.  This  includes  the  effects  of  multi-path,  Doppler 
shifts,  and  fading.  Research  could  be  done  on  how  these  phenomena  affect  probabil¬ 
ity  of  error  in  the  system  and  cross-correlation  of  received  symbols. 

6.3.5  Optimal  Multiple  Access  Codes.  This  thesis  discussed  how  a 
limited  number  of  orthogonal  users  could  be  generated  for  a  given  magnitude  spec¬ 
trum.  Further  research  could  be  done  on  extending  the  number  of  orthogonal  users. 
Another  possibility  is  generating  additional  users  with  symbols  that  have  low  cross 
correlation  with  symbols  of  the  other  network  users;  the  lower  the  cross  correlation 
between  symbols  of  network  users,  the  lower  the  multiple  access  interference. 
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Appendix  A.  Matlab  Simulation  Code 

Matlab  Functions 

function  [ang_vec] =randangl (points , sections) 

°/0This  function  generates  the  random  angles  used  by  the  TDCS. 
/(Variables 

/(points:  the  number  of  points  in  a  symbol.  "points"  needs  to  be  a 
/(power  of  2 . 

/(sections :  the  number  of  different  equally  spaced  angles  used 
°/0ang_vec:  A  vector  of  phases  for  a  basis  function 

if  round(log2(points) )-log2(points) ~=0 

error( ’Number  of  points  needs  to  be  a  power  of  2’); 

end 

if  r ound(log2( sect i ons ) )-log2( sect ions) ~=0 

error ( ’Number  of  sections  needs  to  be  a  power  of  2’); 

end 

if  points<4 

error ( ’Number  of  points  needs  to  be  greater  >=4 ’ ) ; 

end 

°/0random  phases  are  generated  for  the  first  half  of  the  frequency 

“/(components 

halfpoints=points/2-l ; 

A=rand(l , halfpoints) ; 

B=A*sections ; 

C=floor(B) ; 

D=C/ sections*2*pi ; 

/(second  half  of  the  frequency  components  generated  using  the  fliplr 
/(command  and  the 
/(first  half  phases 
ang_vec=[0,D,0,fliplr(-D)]  ; 

function  [basisf ] =basisgenpsk (envelope , angles ,numcodes) 

°/0This  function  generates  the  basis  functions  used  by  the  TDCS. 

/(Basis  functions  are  created  in  the  frequency  domain  and  transformed 
°/0into  the  time  domain. 

/(Variables 
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°/0envelope:  Magnitude  of  the  frequency  components  of  the  basis 
°/0f  unctions 

Wangles :  Phase  of  the  the  frequency  components  of  the  initial  basis 
°/0f  unction 

°/0numcodes:  Number  of  basis  functions 
°/0basisf  :  output  basis  functions 

°/0Note  that  the  magnitude  of  the  frequency  components  for  the  "initial 
°/0basis  function"  are  found  in  the  variable  "envelope"  and  the  phase 
°/0components  in  the  variable  "angles" 

[h,w] =size (envelope) ; 

°/0This  for  loop  adjusts  the  phase  of  initial  basis  function  to 
°/0generate  the  phases  of  the  other  basis  functions 
for  m=l :numcodes 

“/(Magnitude  of  phase  adjustment 
codeangles=ones (1 , w/2-1) * (m-l)/numcodes*2*pi ; 

“/.Vector  of  phase  adjustments 

codeanglet= [0 , codeangles , 0 , f liplr (-codeangles) ] ; 

“/.Phase  adjustment  added  to  phase  of  initial  basis  function 
codeanglev=angles+codeanglet ; 

“/.Vector  is  transformed  to  create  time  domain  vector  of  basis 
“/function 

basisf (m, : )=real (if ft (envelope . *exp( j *codeanglev) ) ) ; 

end 


function  [noiseenv] =env2 (points , signal ,noisef loor) 

“/This  function  generates  a  given  noise  floor  level  based 
°/0on  the  number  of  points  in  a  symbol.  The  length  of  the 
°/0noise  vector  matches  the  length  of  the  signal. 

[h,w] =size(signal) ; 

noiseenv=randn(l , w) * (noisef loor/points) ~ . 5 ; 


function  [grayoutput] =graycoder (bits) 

“/This  function  generates  a  gray  code  of  length  "bits" 
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"/(variable 

°/0bits=length  of  gray  code  word 

°/0grayoutput=matrix  where  each  row  is  a  dif f erentgray  code  word. 
%the  number  of  rows  is  equal  to  2~bits 

A= [0 ; 1]  ; 

"/(Generates  the  case  for  bits=l 
if  bits==l 

grayoutput=A; 

"/(Generates  the  gray  code  for  bits>l 
else 

for  q=l:bits-l 

[c,d]=size(A) ; 

E= [A; f lipud(A)] ; 

F= [zeros (c , 1) ; ones (c , 1)] ; 
grayoutput= [F,E] ; 

A=grayoutput ; 

end 

end 

function  [databits] =datagen(numdatabits) 

"/(Generates  random  vector  of  l’s  and  0’s  of  length  "numdatabits" . 
databits=round(rand(l ,numdatabits) ) ; 

function  [transmit , f inalsymbols] =modulatorb (data, basis) 

°/0This  function  modulates  a  vector  of  ones  and  zeros 
fusing  the  basis  functions  it  is  given. 

/(Variables 

°/0data:  The  data  vector  of  ones  and  zeros 
°/0basis:  A  matrix  where  each  row  is  a  basis  function 
"/(transmit :  transmitted  signal  with  modulated  symbolos 
%f inalsymbols :  Symbols  for  modulation 

"/determines  the  number  of  basis  functions 
[a,b] =size (basis) ; 

"/(determines  the  number  of  data  bits  transmitted 
[c , d] =size (data) ; 

/(concatenates  zeros  to  the  end  of  the  data  string 
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°/0so  the  system  has  an  integer  number  of  symbols  to  modulate 
e=ceil(d/log2(a) ) ; 
f =zeros (log2 (a) ,e) ; 
f (1 :d)=data(l :d) ; 

g=f  ’ ; 

"/generates  gray  code 
log2(a) ; 

h=graycoder (log2 (a) ) ; 

"/assigns  a  symbol  to  groups  of  databits 
k=zeros (1 , e)  ; 
if  a~=2 

for  q=l:e 

[varl , var2] =min(sum( (abs (ones (a, 1) *g(q, : )-h) ) ’ ) ) ; 
k(q)=var2; 

end 

else 

for  q=l:e 

[varl , var2] =min( ( (abs (ones (a, 1) *g(q, : )-h) ) ’ ) ) ; 
k(q)=var2; 

end 

end 

"/assigns  a  basis  function  based  on  the  symbol 
"/basis  functions  are  then  concatenated  one  after  the  other 
f inal=zeros (1 ,b*e) ; 
for  t=l : e 

final ( ( (t-1) *b+l) : (t*b) )=basis (k(t) , : ) ; 

end 

transmit=f inal ; 
f inalsymbols=k ; 


function  [f inalbits ,f inalsymbols] =demodulatorb (transmit , basis) 

"/This  function  demodulates  the  transmitted  signal  using  the 
"/basis  functions  it  is  given. 

"/Variables 

"/transmit:  The  input  signal  from  the  TDCS  transmitter 
"/basis:  A  matrix  where  each  row  is  a  basis  function 
°/0f inalbits:  Demodulated  bits 
°/0f  inalsymbols  :  Demodulated  symbols 
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"/determines  the  number  of  basis  functions 
[a,b] =size (basis) ; 

"/determines  the  size  of  the  transmitted  signal 
[c , d] =size (transmit) ; 

"/determines  the  number  of  symbols  in  the  transmitted  signal 
numsymbs=ceil (d/b) ; 

QQ=zeros (b ,numsymbs) ; 

QQ ( 1 : numsymbs*b) =transmit ; 

RR=QQ  ’  ; 

if  a==2 

A=RR*basis (1 ,  : )  ’  ; 

F=sign(A) *-l ; 
symbolx=( (F+l) /2+1)  ’  ; 


else 

obl=a/4+l ; 

"/breaks  up  transmitted  signal  into  sections  representing  symbol 
A=RR*basis (1 , : )  ’  ; 

B=RR*basis (obi , : ) ’ ; 


C=ones (a+1 , 1) *mod(atan2(B’ , A’ ) , 2*pi) ; 

D=( (2*pi* [0 : a] /a) ’ ) *ones (1 ,numsymbs) ; 

[E , F] =min (abs (C-D) ) ; 

F=F-1; 

F=mod(F,a) ; 
symbolx=F+l ; 

end 

"/generates  a  gray  code  based  on  the  number  of  basis  functions 
h=graycoder (log2 (a) ) ; 

"/assigns  symbol  to  each  section  of  the  transmitted  signal 
"/  finala=zeros(numsymbs,log2(a)) ; 
for  q=l:numsymbs 

f inala(q, : )=h(symbolx(q) , : ) ; 
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end 


f inalb=f inala ’  ; 

f inal=zeros (1 ,numsymbs*log2 (a) ) ; 

final(l :numsymbs*log2(a))=finalb(l :numsymbs*log2(a)) ; 

'/demodulated  bits 
f inalbits=f inal ; 

'/demodulated  symbols 
f inalsymbols=symbolx ; 
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clear  all 
close  all 


iterations=1563; 

numsymbs=64; 

tic 


Bjam= [0.0000,  0.3034,  0.4142,  0.6362,  1.0888,  2.0215,... 

3.4755,  3.8799,  3.6811,  2.9683,  1.9194,  1.2279,. 
0.8696,  0.6794,  0.5667,  0.4928,  0.4406,  0.4017,. 
0.3716,  0.3477,  0.3284,  0.3125,  0.2994,  0.2884,. 
0.2793,  0.2718,  0.2656,  0.2605,  0.2566,  0.2535,. 
0.2514,  0.2502,  0.0000,  0.2502,  0.2514,  0.2535,. 
0.2566,  0.2605,  0.2656,  0.2718,  0.2793,  0.2884,. 
0.2994,  0.3125,  0.3284,  0.3477,  0.3716,  0.4017,. 
0.4406,  0.4928,  0.5667,  0.6794,  0.8696,  1.2279,. 
1.9194,  2.9683,  3.6811,  3.8799,  3.4755,  2.0215,. 
1.0888,  0.6362,  0.4142,  0.3034]; 


for  p=l : 20 

"/.makes  the  notch 

if  p==3  |  p==4  |  p==7  |  p==8  |  p==ll  |  p==12  |  p==15  I  p==16 
I  p==19  I  p==20 

B=[0, ones (1,31) ,0, ones (1,31)] ; 

B=B* (64/sum(B) )  ~  .  5 ; 
varsum=sum(B . ~2) ; 

B(5: 12)=zeros(l,8) ; 

B(54 : 61)=zeros (1 , 8) ; 
varsum2=sum(B . ~2) ; 

B=B*sqrt (varsum/varsum2) ; 

else 

B=[0, ones (1,31) ,0, ones (1,31)] ; 

B=B* (sum(B) /64) ~ . 5 ; 

end 

for  n=l:4 

ary=2~ (n) 

numdatabits2=n*numsymbs ; 
for  m=l : 23 

EbNodb(m)=m-6; 
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EbNo=10~ (EbNodb (m) / 10) ; 

noisef loor=sum(B . ~2) /2/log2 (ary) /EbNo ; 

for  k=l : iterations 

C=randangl (64 , 16) ; 
basisfunc=basisgenpsk(B,C,ary) ; 
Q=datagen(numdatabits2) ; 

[R,R2] =modulatorc (Q ,basisfunc) ; 

if  p==  2  |  p==4  |  p==6  |  p==8  |  p==10  I  p==12 
I  p==14  I  p==16  I  p==18  I  p==20 
Cjam=randangl (64, 128) ; 
basisfunca=basisgenpsk(Bjam,Cjam, ary) ; 
Qa=zeros (1 ,numdatabits2) ; 

[Ra , R2a] =modulatorc (Qa , basisf unca) ; 

else 

Ra=zeros (1 , 64*numsymbs) ; 

end 

if  p>4 

if  p>12 

D=C (2 : 32) ; 

F=zeros(l,31) ; 

F ( 1 , : ) =D ; 
qxx=l ; 

for  pxx=l:31 

F (pxx) =D (pxx) +2*pi*pxx*qxx/31 ; 

end 

Cb=[0,F,0,f liplr(-F)] ; 


else 

Cb=randangl (64, 16) ; 

end 

basisf uncb=basisgenpsk (B , Cb , ary) ; 
Qb=datagen(numdatabits2) ; 

[Rb , R2b] =modulatorc (Qb , basisf uncb) ; 
if  p==  9  |  p==10  I  p==ll  I  p==12  |  p==17 
p==18  |  p==19  I  p==20 
tau=ceil (rand(l) *64) ; 
if  tau==64 
Rb=Rb ; 
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else 


Rb= [Rb(tau+1 : length (Rb) ) ,Rb(l :tau)]  ; 

end 

else 

Rb=Rb ; 

end 

else 

Rb=zeros (1 , 64*numsymbs) ; 

end 


noisestuf f =env2 (64,R,noisef loor) ; 
S=R+Ra+Rb+noisestuf f ; 

[T,T2] =demodulatord(S ,basisfunc) ; 
per cent error sx (k) =sum (sum (abs (T ( 1 : numdatabits2) -Q 
) ) ) / (numdatabits2) ; 

percentsymbolerrorsx(k)=sum(sum(abs (sign(R2-T2 . . . 
) ) ) ) /length (R2) ; 


end 

percenterrors (m , n , p) =mean (percenterrorsx) ; 

percent symbolerrors (m , n , p) =mean (percentsymbolerrorsx) 

end 

end 

end 

toe 

save  btest21  EbNodb  percenterrors  percentsymbolerrors  iterations 
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°/0Abel  Nunez 
°/„EENG  799 

/(Multiple  Access  effects  on  single  user  by  a  network  of 
°/0orthogonal  users 

clear  all 

close  all 

EbNo=6 ; 

kary=2 

ary=2~kary ; 

numdatabits2=64*kary ; 

points=64;  °/„length  of  symbol 

iterations=15625;  %number  of  iterations  the  program  is  run 
tic 

7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 

/  o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /  0  /o 

pberrorf in=zeros (16 , 32) ; 
pserrorf in=zeros (16 , 32) ; 
counter=0 ; 
for  px=l : 16 
pp=px-l; 

pberrortot=zeros (1 ,32) ; 
pserrortot=zeros (1 ,32) ; 

°/„  jamming 

if  mod(pp,2)==l 

Bjam= [0.0000,  0.3034,  0.4142,  0.6362,  1.0888,  2.0215,... 

3.4755,  3.8799,  3.6811,  2.9683,  1.9194,  1.2279,. 

0.8696,  0.6794,  0.5667,  0.4928,  0.4406,  0.4017,. 

0.3716,  0.3477,  0.3284,  0.3125,  0.2994,  0.2884,. 

0.2793,  0.2718,  0.2656,  0.2605,  0.2566,  0.2535,. 

0.2514,  0.2502,  0.0000,  0.2502,  0.2514,  0.2535,. 

0.2566,  0.2605,  0.2656,  0.2718,  0.2793,  0.2884,. 

0.2994,  0.3125,  0.3284,  0.3477,  0.3716,  0.4017,. 

0.4406,  0.4928,  0.5667,  0.6794,  0.8696,  1.2279,. 

1.9194,  2.9683,  3.6811,  3.8799,  3.4755,  2.0215,. 

1.0888,  0.6362,  0.4142,  0.3034]; 

else 

Bjam=zeros(l,64) ; 

end 

/(notching 
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if  mod(f loor (pp/2) ,2)==1 

B= [0 , ones (1 , point s/2-1) , 0 , ones (1 ,points/2-l)] 
B=B*(sum(B)/64)~ .5; 
varsum=sum(B . ~2) ; 

B(5: 12)=zeros(l,8) ; 

B(54 : 61)=zeros (1 , 8)  ; 
varsum2=sum(B . ~2) ; 

B=B*sqrt (varsum/varsum2) ; 

else 

B= [0 , ones (1 , point s/2-1) , 0 , ones (1 ,points/2-l)] 
B=B* (sum(B) /64) ~ . 5 ; 

end 

noisef loor=sum(B . ~2) /2/log2 (ary) /EbMo ; 
if  px==5  |  px==6 

iterations=1563*10 ; 

else 

iterations=1563; 

end 

for  iteration=l : iterations 


TIT/  TH 

/o  /o  /o  /o  /o  /o  /  0 

“/.Generate  Angles 

V  o/o/ o/o/ 0/0/ 

/ 0  /o  /o  /o  /o  /o  /o 

if  mod(f loor (pp/4) ,2)==1 
%orthogonal  angles 
Co=randangl (64 , 16) ; 

D=Co (2 : 32) ; 

F=zeros(16,31) ; 

F ( 1 , : ) =D ; 

for  q=l:30 

for  p=l : 31 

F (q+1 , p) =D (p) +2*pi*p*q/31 ; 

end 

end 

G=  [zeros (31,1) , F, zeros (31,1) , f liplr (-F)] ; 
G= [G; zeros (1 , 64)] ; 
for  q2=32:32 

G (q2 , : ) =randangl (64,16); 

end 

C=G; 

else 
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"/regular  angles 
G= [zeros (32, 64)] ; 
for  q2=l:32 

G (q2 , : ) =randangl (64,16); 

end 

C=G; 

end 

for  kk=l:32 

basisfunc( : , : ,kk)=basisgenpsk(B,C(kk, : ) ,ary) ; 

end 

pberrorsx=zeros (1,32) ; 

Q=datagen(numdatabits2) ; 

[R,R2] =modulatorc (Q , basisf unc ( : , : , 1) ) ; 

C j  am=randangl (64,128) ; 

basisf unca=basisgenpsk(Bj am , C j  am, ary) ; 

Qa=zeros (1 ,numdatabits2) ; 

[Ra,R2a] =modulatorc(Qa, basisf unca) ; 

Interf erence=env2 (points ,R,noisef loor) ; 

S=Interf erence+R+Ra ; 

for  users=l:32 

[T,T2] =demodulatord(S , basisf unc ( : , : , 1) ) ; 
pberrorsx(users)=sum(sum(abs (T(l : numdatabits2) -Q) ) )  . . 
/ (numdatabits2) ; 

pserrorsx(users)=sum(sum(abs(sign(R2-T2) ) )) /length (R2) 
if  users~=32 

Qa=datagen(numdatabits2) ; 

[Rb,R2b] =modulat or c(Qa, basisf unc ( : , : ,users+l) ) ; 

if  mod(f loor (pp/8) ,2)==1 
tau=ceil(rand(l)*64) ; 
if  tau==64 
Rb=Rb ; 

else 

Rb= [Rb(tau+1 : length (Rb) ) ,Rb(l : tau)] ; 

end 
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else 


Rb=Rb ; 

end 


S=S+Rb ; 

end 


end 

pberrortot=pberrortot+pberrorsx; 
pserrortot=pserrortot+pserrorsx; 
if  mod(iteration, 100)==0 
count er=counter+l 

end 

end 

pberrorf in(px, : )=pberrortot/ iterations ; 
pserrorf in(px, : )=pserrortot/ iterations ; 

end 

toe 

save  bmanet21  pberrorf in  pserrorf in  EbNo  iterations 
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°/0Abel  Nunez 
°/„EENG  799 

°/0Multiple  Access  effects  on  single  user  by  a  network  of 
°/0orthogonal  users 


clear  all 
close  all 

points=64;  °/„length  of  symbol 

iterations=1563;  °/„number  of  iterations  the  program  is  run 

B=  [0 , ones (1 , point s/2-1) , 0 , ones (1 ,points/2-l)] ; 

/(frequency  envelope.  DC  and  Nyquist  frequency  are  nulled  out 

counter=0 ; 

EbNo=6 ; 
kary=2 
ary=2~kary ; 

noisef loor=sum(B . ~2)/2/log2 (ary)/EbNo ; 
numdatabits2=64*kary ; 

tic 

pberrortot=zeros(l,31) ; 
pserrortot=zeros(l,31) ; 

for  iteration=l : iterations 
pberrorsx=zeros (1 ,  31) ; 
pberrorsx=zeros (1 ,  31) ; 
for  mm=l : 5 

Cx=randangl (64, 16) ; 

basisfunc( : , : ,mm)=basisgenpsk(B,Cx,ary) ; 
Q=datagen(numdatabits2) ; 
if  mm==l 
Qx=Q; 

Cc=Cx ; 

end 

[R(mm, : ) ,R2 (mm, : )] =modulatorc (Q ,basisfunc ( : , : ,mm) ) ; 

end 

Rtot=sum(R) ; 
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for  notcher=l:31 

BN=[0,ones(l,32-notcher) , zeros (1 ,notcher-l) ,0, . .  . 

zeros(l ,notcher-l) ,ones(l ,32-notcher)] ; 

BN2=BN.*B; 

basisf unc2=basisgenpsk(BN2 , Cc , ary) ; 
noisef loor=sum(BN2 . "2) /2/log2 (ary) /EbNo ; 

Interf erence=env2 (points ,Rtot , noisef loor) ; 

S=Interf erence+Rtot ; 

[T,T2] =demodulatord(S ,basisfunc2) ; 

pberrorsx(notcher)=sum(sum(abs(T(l :numdatabits2)-Qx) ) ) . . . 
/ (numdatabits2) ; 

pserrorsx(notcher)=sum(sum(abs(sign(R2(l , :)-T2)))) . . . 
/length(R2 (1 , :)) ; 


end 

pberrortot=pberrortot+pberrorsx; 

pserrortot=pserrortot+pserrorsx; 

end 

pberrorf in=pberrortot/ iterations ; 
pserrorf in=pserrortot/ iterations ; 

toe 

save  bnotchl  pberrorf in  pserrorf in  EbNo  iterations 
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°/0Abel  Nunez 
°/„EENG  799 

°/0Multiple  Access  effects  on  single  user  by  a  network  of 
°/0orthogonal  users 


clear  all 
close  all 

points=64;  °/„length  of  symbol 

iterations=1563;  °/„number  of  iterations  the  program  is  run 

B=  [0 , ones (1 , point s/2-1) , 0 , ones (1 ,points/2-l)] ; 

/(frequency  envelope.  DC  and  Nyquist  frequency  are  nulled  out 

counter=0 ; 


EbNo=6 ; 
kary=2 
ary=2~kary ; 

noisef loor=sum(B . ~2) /2/log2 (ary)/EbNo ; 
numdatabits2=64*kary ; 

tic 

pberrortot=zeros(l,31) ; 
pserrortot=zeros(l,31) ; 

for  iteration=l : iterations 

Co=randangl (64 , 16) ; 

D=Co (2 : 32) ; 

F=zeros(30,31) ; 

F ( 1 , : ) =D ; 

for  q=l:30 

for  p=l : 31 

F(q+1 ,p)=D(p)+2*pi*p*q/31 ; 

end 

end 

G=  [zeros (31,1) , F, zeros (31,1) , f liplr (-F)] ; 
G= [G; zeros (1 , 64)] ; 
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for  kk=l : 30 

basisf unc ( : , : , kk) =basisgenpsk (B , G (kk , : ) , ary) ; 

end 


for  mm=l : 5 

Q=datagen(numdatabits2) ; 
if  mm==l 
Qx=Q; 

end 

[R(mm, : ) ,R2 (mm, : )] =modulatorc (Q , basisf unc ( : , : ,mm) ) ; 

end 

Rtot=sum(R) ; 

pberrorsx=zeros (1 , 31) ; 
pberrorsx=zeros (1 , 31) ; 
for  notcher=l : 31 

BN=[0,ones(l,32-notcher) , zeros (1 ,notcher-l) ,0, . . . 

zeros(l ,notcher-l) ,ones(l ,32-notcher)] ; 

BN2=BN.*B; 

basisfunc2=basisgenpsk(BN2,G(l , : ) ,ary) ; 
noisef loor=sum(BN2 . "2) /2/log2 (ary) /EbNo ; 

Interf erence=env2 (points ,Rtot , noisef loor) ; 

S=Interf erence+Rtot ; 

[T,T2] =demodulatord(S ,basisfunc2) ; 

pberrorsx(notcher)=sum(sum(abs(T(l :numdatabits2)-Qx) ) ) . . . 
/ (numdatabits2) ; 

pserrorsx(notcher)=sum(sum(abs(sign(R2(l , :)-T2)))) . . . 
/length (R2 (1 , :)) ; 


end 

pberrortot=pberrortot+pberrorsx; 

pserrortot=pserrortot+pserrorsx; 

end 

pberrorf in=pberrortot/ iterations ; 
pserrorf in=pserrortot/ iterations ; 

toe 

save  borthonotchl  pberrorf in  pserrorf in  EbNo  iterations 
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Appendix  B.  Cross-Correlation  Between  PSK  Symbols 

n=time  index 

N=number  of  samples  that  make  up  a  symbol. 

P= number  of  sinusoids  which  make  up  a  symbol.  P  <  y 

0cU)=phase  modulation  of  cth  symbol  for  vth  user.  S$  —  |o,  ...  ,  ('Q~q27T  | 

(fin"1  =phase  coding,  superscript  indicates  a  set  of  phase  coding  for  vth  user 

Symbol  1  is  defined  with  the  following  equation.  Superscript  indicates  user. 
Subscript  indicates  symbol. 

4v)  (n)  =  Ap  cos  ( 2rrP ^  +  9c})  (B-1) 

p= i 

Symbol  2  is  defined  with  the  following  equation.  Superscript  indicates  used. 
Subscript  indicates  symbol. 

sdV)  ( n )  =  Bp  cos  (27rP^  +  4U,)  +  9P)  (B-2) 

p=  i 

B.l  Cross  Correlation  Between  Symbol  1  and  Symbol  2 


N- 1 

R\ P  =  E  Sf  (n)  S<”>  («) 

n=0 


(B.3) 


AT-i  P  P 

COS  f2^  +  ^  +  9c])  COS  (2?r^  +  ^  +  9> 


N 
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(w) 

d 
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P  P  N-l 


N 2 


EE  EE  E] apbv cos (27Tpjj  +  4u)  +  0cv)) c°s (2n^^  +  ^  (B-5) 


p=  1  g=l  n=0 


'AT 


P  P  N-l 

E  E aa cos {2npjj + 4W) + A}) cos (27r^ + ^ + °{d] 

p=  1  q=l,q^p  n= 0 
4  p  iV"1 


+  ]v2  EE  EE  a  a- cos  (27rr^  +  ^  +  9^) cos  (27rr^  +  ^  +  9*w)) 


r=l  72=0 
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P  P 
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P  P 


TV— 1 
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EE  EE  AP59EE2cos(27r^_g)^+^)_^)+^)_^w)) 


p=l  q=l,q^p  n= 0 

4  ^  ^  1 


Y  ArBr  Y  2  cos  (4^  +  A4  +  <A}  +  A4  + 


r=l  72=0 

4  A  ^  i 

EE  A-^r  EE 

r=l  ?2=0 


—2  Y  nb,  Y  72  cos  (4”’  -  4”’  +  -  4f ’) 
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A  A  9  cos 
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,M  _  +  gM  _  #(”) 
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Ea 
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B.2  Cross  Correlation  Between  Symbols  of  a  Single  user 

v  =  w 


R 


(v,v) 

cd 


(B.10) 


R 


(v,v) 

cd 


2 

N 


Al  cos 


p= i 


(B.ll) 


B.3  Autocorrelation  of  symbol 


td(v,v) 

1  Lcc 


2 

N 


p 

Al  cos (°) 

p= i 


(B.  12) 


P=  1 


(B-13) 
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Mean  of  cross-correlation  between  symbols  of  two  different  users.  Each  symbol 
is  assumed  to  be  equiprobable,  v  ^  w,  and  users  v  and  w  are  synchronized  (symbol 
boundaries  coincide). 


*Sr)  =  ^  E  aprp  cos  ($°  -  ^ + w  -  dP) 

p=i 


(B-14) 


Mean 


■p  \  d{v,w) 

^0  I  Kcd 


=  E*  \  ^  E  cos  ($°  -  4")  +  °iv)  - 

p= i 


(B.15) 


P 


-  £  ApJBpE,  {cos  ($>>  -  -  0?° 

p= i 


)}  (B.16) 


=  4  E  4^®*  {cos  -  ^w>) «*>  (4’1  -  4”’) } 

p=i 

-4  E  44E,  {sin  (4”  -  4”’)  si»  «’  -  4”)) }  (B-17) 

P=1 

=  Jr  E  44E,  {cos  -  4”1)  }  cos  «>  -  <"') 

P=1 

-4  E  ArB- E»  {sin  (»?”  -  4”’)  }  sin  «>  -  4”’)  (B.18) 

P=1 


B-4 


=  X  ApRpEo  {cos  (ecj)) cos  ( 6dw) ) }  c°s  W  -  4w)) 

p= i 

{sin  (^}) sin  (0aw)) } cos  Ww)  _  4w)) 

p=i 

E  vy=#  {™  (0<”>)  cos  («<"'>)  }  sin  «>  -  <”>) 

P=1 

X  AAEe  {cos  (&cv))  sin  (&dW))  }  sin  {(pp]  -  <PA))  (B.19) 

p=  i 


|  E  4cBpE»  {cos  (OW) }  E„  {cos  («<’">)  }  cos  «■>  -  $*>) 

p=l 

+  ^  X  Aspe*  (sin  (#cv))  }  E0  {sin  (X'°)  }  COS  (<X  -  4w)) 

p=i 

-|x  aae<?  (sin  (A°) } E^  {c°s  (0dw)) } sin  (X  -  4U°) 

p=i 
2  P 

-J7  E  4pBpE„  {cos  («<’>) }  E„  {sin  («<”>)  }  sin  (*W  -  $">)  (B.20) 

P=1 

E«{Pir’}=0  (B.21) 
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Variance  of 
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jy: 2  X/  ApBpAqBq 

P=  1  9=1 

xE0  {cos  (4")  -  +  0^  -  +  0W  -  ^})  } 

P  P 

+  jy2  ApBpAqBq 

P=  1  9=1 

xEo  {cos  (<^  -  <j>W  +  0<«)  _  0^  -  $)  +  -  0(«)  +  0^)  }  (B.26) 


jy2  ApBpAqBq 

P=  1  9=1 

xEo  {cos  -  0^  +  0<w>  -  0M  +  ey  -  0H  +  0W  -  0<w))  } 

P  P 

+  jy2  \BpAqBq 

P=  1  9=1 

xEd  |cos  0<”>  -  4”>  -  ./.<*>  +  $">  +  0W  -  9'”1  -  «<•>  +  «<”>)  }  (B.27) 
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]\[2  ^22  22  ApBpAqBq 
P=  1  9=1 

xE„  j  cos  ($>>  -  4”>  +  0<’>  -  </.<«>)  cos  (OW  -  »<•">  +  9'”1  -  «?’)  } 

P  P 

_  ]y2  ApBpAqBq 

P=  1  9=1 

xEoj  sin  (*W  -  4”>  +  0W  -  *<”>)  sin  (o<”>  -  0<”>  +  9<«>  -  9*"1)  | 

P  P 

+  ]y2  ApBpAqBq 

P=  1  9=1 

XE„  j  cos  (0W  -  4”>  -  *<«>  +  <#.<”>)  cos  (««”>  -  9<”>  -  9<”>  +  9<”>)  | 
p  p 

-  7V2  ApBpAqBq 

P= 1  9=1 

xE„|  sin  «>  -  <”>  -  4«>  +  4»>)  Sill  (»<«>  -  0<”>  -  #W  +  #<•">)  |  (B.28) 
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p  p 

=  ^EE  ArBpA,B,  cos  «>  -  4”>  +  <$)  -  *<”>) 

P=1  9=1 

xE,  {cos  (of  -  «f  +  Of  -  «f)  } 

P  P 

-&ZE  ApBpAqBq  sin  (0j0}  -  0jf)  +  0$u)  -  0^}) 

P=  1  9=1 

xE,  {sin  (of  -  of  +  Of  -  of )  } 

o  p  p 

+^EE  WW  cos  «)  -  4”>  -  4«>  +  <”>) 

P=1  9=1 

xe,  {cOS(of -of -of +  of)} 

P  P 

-0*EE  ApBpAqBq  sin  (0^  -  ^  -  0jf)  +  0£u)) 

P=1  9=1 

xE,  {sin  (of  -  Of  -  Of  +  Of)  }  (B.29) 
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-fEE  ApBpAqBq  sin  (0^  -  0^  +  0^  -  0^)  Ee  {sin  (20^  -  20<w))  } 


p=i  9=1 

p  p 


t0EE  ApBpAqBq  cos  f f  -  -  0W  +  0f )  E„  {cos  (0)} 


AT2 

2 

'7V2 


p=l  9=1 
p  p 


EE  ApBpAqBq  sin  (0^  -  0j°>  -  0^  +  0^)  Ee  {sin  (0)}  (B.30) 


p=i  9=1 
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2 

jY 2  ApBpAqBq  COS 

P=1  9=1 

P  P 

-  ]V2  ApBpAqBq 

P=  1  9=1 
P  P 

+  jy2  ApBpAqBq 

P=  1  9=1 


«’  -  4”’  +  4”’  -  <”’)  E»  {cos  (2 

sin  -  cJj,"1  +  0<">  -  0*"1)  E«  {sill 
cos  (</.<”>  -  $*>  -  <'  +  </.(«>) 


(B.31) 


P  P 

=  ^EE  ApBpA,B,  cos  (</.<«>  -  $">  +  -  </.<”>) 

P=1  9=1 

xE0  |  cos  (2div))  cos  (29(^  +  sin  ( 29(cv) )  sin  (26^  } 

P  P 

-^EE  ApBrA,B,  sin  ($”>  -  +  0<”>  -  <”>) 

P=1  9=1 

xE„  | sin  (2 9^)  cos  (2^)  -  cos  (20^)  sin  (29^  } 

P  P 

+^EE  ApBpA,B,  cos  «”>  -  4«>  -  +  </><”>)  (B.32) 

P=1  9=1 
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p  p 

=  ^£E  ArBpA,B,  cos  «>  -  4”>  +  4”>  -  4”>) 

P=1  9=1 

xE0  {cos  (2^))}E,{cos  (2^)} 

P  P 

+]0££  ApBpAqBq  cos  (4”)  -  4”>  +  0<”)  -  4”>) 

P=  1  9=1 

xE0  {sin  (20^)  }  Efl  {sin  (20<w))  } 

P  P 

-^££  ApBpAqBq  Sill  (4”>  -  4”>  +  4”>  -  4”)) 

P=1  9=1 

xE,{sin(2^)}E,{cos  (2^)} 

P  P 

+0s££  ApBpA,Bq  Sill  (4”)  -  4”>  +  4”>  -  4”)) 

P=1  9=1 

xEe  {cos  (2»W)}E,{sin  (2#'"1)} 

P  P 

q0££  ApBpAqBq  cos  (0^  -  0^}  -  0^  +  0^)  (B.33) 

p=i  9=1 


var 


{^r’ }  =  Jj  £  £  ApBpAqBq  cos  (4«)  -  4»>  +  4«)  -  4”») 

p=l  9=1 

xE,{cos  (2^))}E,|cos  (20^)} 

P  P 
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(B.34) 


Variance  of  R^’p'> 


in  terms  of  0  for  BPSK.  This  equation  describes  the  variance 


with  no  restrictions  on  0. 
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p  p 


var  [B{^dw) }  =  ^2  Y1  Y1  ApBp\b<i  cos  W  ~  4W)  +  “  ^}) 


P=1  9=1 


P  P 


+^EE  ApBpAqBq  cos  (<^  -  4U,)  -  <f)W  +  <^>)  (B.35) 


P=1  9=1 


P  P 


var  {  ^r} }  =  ]^2  Z)  S  ApBpAiB‘ 


))cos(^-^)  (B.36) 


9  COs(0£°  -  0j"})  COs(0g1')  -  0jo)'' 


P=1  9=1 


Variance  of 


in  terms  of  <f>  for  M 


2k  PSK  where  k  >  1.  This  equation 


describes  the  variance  with  no  restrictions  on  0. 


vaWdT’V^EX> 

P=1  9=1 


2 

n2 


pBpAqBq  COS 


-  V*  -  V  +  V)  <B-37) 


The  reason  (B.36)  and  (B.37)  differ  is  the  expected  value  of  cos(20)  is  1  for 
BPSK  and  0  for  MPSK,  where  M  =  2k ,  k  >  1. 
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B.4  Expected  Value  of  the  Variance  for  the  BPSK  Case  with  Random 
Data  Phase  Modulation 

Assuming  cf)  is  uniformally  distributed,  U[0,27t),  the  expected  value  for  the 
BPSK  case,  (5.36)1,  becomes 


E A  (5.36)  \  =E 


p 


^££  A,BvA,B ,  cos«>  -  4“>)  cos(^>  -  <”>) 

P=  1  9=1 


(B.38) 


p  p 

£  £  ApBpAqBq  COS (0^}  -  0j"})  COs(0^} 
P=1  9=1, 


+E^ 


(B.39) 


y'  .  i  i 

=  MA££  A„BrA„Bt cos(4”>)  cos(4”))  cos«>  -  0<»>) 

^  p=l  <J=1 


(B.40) 
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4EE  ApBpAqBqEj  cos(0^)|  cos(^)  cos(^  -  0j>}) 

p=  1  q= 1  ^  J 

a  p  p  r  i 

-4EE  4,£p.4g£gE  J  sin(0<w>)  l  sin^M)  cos(0^  -  0j>}) 

p=l  q=  1  ^  J 

+^I  E  ^  Q  +  ^E*{  »s(2<’)  -  H”1) }) 

(B.41) 

=  ^E  A2-B ’  G + G*{  cos(2^)  “  2^>)}) 

=  ^E^G) 

p= i  v  7 

Ei(B.36)}  =  -|j^^Br2  (B-42) 

^  7  p=l 
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B.5  Expected  Value  of  the  Variance  for  the  MPSK  Case  with  Random 
Data  Phase  Modulation 

Assuming  cf)  is  uniformally  distributed,  U[0,27t),  the  expected  value  for  the 
MPSK  case.  (£>.37)1,  becomes 
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iV2 


P  P 
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P  P 
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P=1  9=1-97 
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V2 


(B.44) 
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p  p 
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2 

iV2 
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(B.45) 
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(B.46) 


W)}  =  ^X>tB 

t  )  r=1 


2 

r 


B.6  Expected  Value  of  the  Variance  for  the  BPSK  and  MPSK  Case  with 
Random  Multiple  Access  Phase  Coding 

if  <p  is  uniformally  distributed,  U[0,27t),  the  expected  value  for  the  BPSK  case 
is  the  same  as  the  MPSK  case 


2 

iV2 


E  ArBr 


(B.47) 
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B.7  Orthogonal  Users  for  PSK  Symbols 

Symbol  c  is  defined  with  the  following  equation.  Superscript  indicates  user. 
Subscript  indicates  symbol. 
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p 

'Y,  Ap  cos  (^2ttp 
p= i 


n 
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v  =  0, ...,  G  —  1 


(B.48) 


=  (Pp  +  Ys  gn(A)^“ 

<2=1 

P 

G  =  Y  sgn(A-) 

r= 1 


/ 

sgn(t)  =  < 


1 

0 

-1 


t  >  0 


t  =  0 


t  <  0 


G=number  of  users.  ^  is  a  positive  integer. 

Symbol  d  is  defined  with  the  following  equation.  Superscript  indicates  user. 
Subscript  indicates  symbol. 
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Cross  correlation  between  Symbol  c  and  Symbol  d 


dr,=£5t’(^R) 


(B.50) 


rT  =  |f  x>? cos  (<’  -  <” + ^ 


(B.51) 
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A)  If  v=w 


P=i 


B)  If  v  7^  w  and  Ap  can  only  take  on  values  of  A  or  0 


(B.53) 
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